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ABSTRACT 


Recent studies have shown that various inaccuracies occur when 
a single-degree-of-freedom gyroscope is subjected to an angular vibra- 
tory environment. Under such conditions the performance of the gyro- 
scope has been analyzed and the nature of the inaccuracies investigated. 
Errors due to “coning", "dynamic unbalance”, and "sculling™ have been 
included in the investigation. 


Values of the inaccuracies representative of those expected from 
tests using the precision angular vibrator developed by the Instrumentation 
Laboratory of the Massachusetts Institute of Technology are presented. 


By reference to a dynamic model the equation of motion was ob- 
tained. A suitable solution to the equation was determined in series 
form for the case of sinusoidal inputs simultaneously applied about two 
axes, No restrictions are placed on the frequency of the applied mo- 
tions. The terms in the series solution reveal the nature of the angular 
displacement of the float with respect to inertial space about the output 
axis for two cases: (1) zero elastic restraint, (2) finite elastic restraint. 


It is strongly recommended that correlation of actual experimen- 
tal results with the developed theory be obtained prior to further 
theoretical considerations. 


Thesis Supervisor: Dr. Walter Wrigley 


Title: Professor of Instrumentation and 
Astronautics 
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OBJECT 


The object of this thesis is to analyze the performance 
of a floated single-degree-of-freedom integrating gyroscope 
subjected to angular vibrations in order to determine the 


inaccuracies that may result from such an environment. 





CHAPTER 1 


INTRODUCTION 


The age of the golden sixties is upon us. Ushered in with this 
decade is an era of accomplishment which has raised the importance 
of the gyroscope to an all-time high. This era of great advancement 
includes those fields of Science and Engineering which are brought to 
the fore when one thinks of space. Because of space, these fields 
have been searched diligently for methods of navigation and control 
and from this searching the gyroscope has emerged, holding a posi- 


tion of extreme prominence. 


With this relatively new stature, the gyroscope has naturally 
been investigated in efforts to improve its performance as well. as 
to predict what may be expected in the future. The floated single- 
degree-of-freedom integrating gyroscope, an important element in 
most space navigation and guidance systems, is the subject of this 
analysis. The SDFI gyroscope has been interrogated from the 
standpoint of being subjected to varying vibratory inputs which may 
introduce error. These vibrations may be of very low frequency 
due to the bending modes of long missile bodies or they may be of 
high frequencies associated with aerodynamic flutters and buzzing. 


Through the theoretical results, it is hoped that undesirable effects 





caused by such vibratory motions will be recognized and later minimized 


when actual verifying data is available. 


The floated single-degree-of-freedom integrating gyroscope is 
an instrument capable of sensing angular velocities to a high degree 
of accuracy. When utilized in inertial guidance systems, the gyro- 
scope must have considerable stability and extreme sensitivity. The 
vibratory environment in which the gyro functions may introduce un- 
wanted angular outputs. This degradation of desired performance is 
the result of both linear and angular oscillations. An analysis of the 
undesired effects is made in order to determine the nature and magni- 


tudes of the errors involved. 


This investigation was conducted with the view that the Massachu- 
setts Institute of Technology Instrumentation Laboratory will conduct 
tests on a precision angular vibrator. A previous investigation has 
disclosed what may be expected if two equal frequencies are simul- 
taneously applied to the output and spin reference axes of a typical 
single-degree-of-freedom gyroscope. The present analysis is concerned 
with any two frequencies simultaneously applied to the output and spin 


reference axes of a SDFI gyroscope. 


The existing equation of motion is based on a simplified dynamic 
model which accounts for damping, compliance, and inertia. As no 
closed form solution to the equation is known, a series form was used 
which gives a very clear representation of the gyroscopic drift terms. 
The series solution appears easy to handle and analyze in view of the 
theory of superposition. The terms seem reasonable in magnitude and 
should be easily verified by the Massachusetts Institute of Technology 


angular vibrator. 





CHAPTER 2 


THE FLOATED SINGLE-DEGREE -OF-FREEDOM 
INTEGRATING GYROSCOPE 


2.1 Design Features 


The types and locations of components within the ideal floated 
single-degree-of-freedom integrating gyroscope unit is pictorially 
shown in Fig. 2-1. These components are the gyroscopic element, 
gimbal, viscous damper, microsyn torque generator, microsyn signal 


generator, gimbal bearings, and the case. 


The gyroscopic element and gimbal are hermetically sealed in 
a cylindrical shell which is separated from the outer case by a high 
density fluid. The case is also hermetically sealed. The two shells 
form a thin-walled cylinder that performs two vital functions when 
completely filled with a high density fluid: 

1. it provides a viscous damping torque that opposes the 


precessional torque generated by the gyro element, 


2. it provides a buoyant force that neutralizes the weight 


of the gyro element and its supporting structure. 


The clearance between the gimbal float and the case plus the fluid 
viscosity determines the damping of the viscous damper; the fluid 


density and the amount of displaced fluid determine the buoyant force. 





Power is conducted to the hysteresis type alternating-current 
motor by flex leads. The entire float assembly, or torque summing 
member, is initially aligned with the case by jewel bearings. Final 
alignment is by magnetic suspension. The suspension (elastic re- 
straint to float radial motions) is provided by the microsyns in con- 
junction with their primary functions of angle measurement and torque 
generation. The magnetic forces necessary for suspension originate 
with the change in excitation currents of properly tuned microsyn 
coils as a result of the radial displacement of the microsyn rotor from 
its centered position. Figure 2-2 is a sectioned view of a typical 
floated single-degree-of-freedom integrating gyroscope. A more de- 
tailed discussion of the single-degree-of-freedom integrating gyroscope 


may be found in Refs. 1 and 2. 
2.2 Operation 


An ideal single-degree-of-freedom integrating gyroscope receives 
as inputs, the angular velocity of its case with respect to inertial space 
about the input axis and electrical current input signals to the torque gen- 
erator. The unit generates an output voltage proportional to the inputs. 
In open loop operation, the gyro output signal would represent the angu- 
lar deviation of the gyro case from a reference orientation. The refer- 
ence is established from initial conditions and by command signals to 


the torque generator. 


In operation, the gyro unit is housed in a thermally controlled 
environment that maintains the temperature to within a few degrees 
Fahrenheit. This control is important because temperature affects the 


density and viscosity of the fluid. 
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OUTPUT AXIS 






GYROSCOPIC aw 
ELEMENT a (OA) 
(ge) ae FIXED TO CASE 
< PARALLEL TO 
GIMBAL GIMBAL BEARING 
SIGNAL (gim) 
GENERATOR 
AND MAGNETIC 
SUSPENSION UNIT 
(ea) ins) SPIN REFERENCE AXIS 
\ (SRA) 
FIXED TO CASE AT RIGHT 
GIMBAL SN ne cone ee 
BEARING SX mas SPIN AXIS OF THE GYRO- 
(gb) ~ SCOPIC ELEMENT WHEN THE 
Ss a OUTPUT SIGNAL FROM THE 
A SIGNAL GENERATOR IS AT 
| ITS NULL LEVEL 
by VISCOUS 
i ve SPIN AXIS GIMBAL ANGLE 
ca 
{ FIXED TO GIMBAL, A : 
INPUT AXIS PARALLEL TO AXIS a oan 
(1A) ABOUT WHICH GYRO er (SRA) 
FIXED TO CASE AT Rue oents 
RIGHT ANGLES TO 
THE OUTPUT AXIS aN 
AND THE SPIN REF- SN REFERENCE DIRECTION — 
ERENCE AXIS deveaewer Coe ees 
DIRECTION 
ae ~ A\(ref)-(empd)] 
~ wa 
NOTES: 1: THE POSITIVE SENSES SHOWN BY THE ARROWS ARE CHOSEN SO THAT (IA), (OA), / 
AND (SRA) FORM A RIGHT-HANDED SYSTEM. Satie 
2. THE SYMBOL. Al (ref)-(empd)] REPRESENTS THE ANGLE A MEASURED FROM THE DIRECTION 
REFERENCE DIRECTION (DENOTED BY (ref) IN THE SUBSCRIPT) TO THE COM- (cmpd) 


PARED DIRECTION (DENOTED BY (cmpd) IN’ THE SUBSCRIPT).* 


DEFINITIONS: 
CASE, (ca) or (c) GYRO UNIT, (gu) 
THE STRUCTURE THAT GIVES:‘SUPPORT FOR THE THE ENTITY MADE UP OF THE COMPONENTS REPRESENTED 
INTERNAL WORKING PARTS OF THE GYRO UNIT, IN THIS DIAGRAM AND ALL THE ADDITIONAL PARTS NECES- 
ENCLOSES THESE PARTS, AND CARRIES PROVI- SARY FOR A SINGLE PACKAGE TO CARRY OUT THE FUNC 
SIONS FOR EXTERNAL CONNECTIONS OF ALL KINDS. TIONS OF A GYRO UNIT. 
TORQUE GENERATOR AND MAGNETIC SUSPENSION SIGNAL GENERATOR AND MAGNETIC SUSPENSION UNIT, 
UNIT, (tg) + (ms) (sg) + (ms) 
COMPONENT FOR (a) RECEIVING INPUT SIGNALS AND COMPONENT FOR (a) RECEIVING THE ANGLE OF THE SPIN 
PRODUCING A CORRESPONDING OUTPUT TORQUE AP- AXIS WITH RESPECT TO THE CASE AS ITS INPUT AND PRO- 
PLIED TO THE GIMBAL ABOUT THE OUTPUT AXIS AND DUCING A CORRESPONDING SIGNAL THAT SERVES AS THE 
(b) ALSO PROVIDING MAGNETIC FORCE ACTING TO OUTPUT SIGNAL FROM THE GYRO UNIT AND (b) ALSO PRO- 
CENTRALIZE THE PIVOT SHAFT IN THE BEARING. VIDING MAGNETIC FORCE ACTING TO CENTRALIZE THE 
DAMPER, (dmp) or (d) PIVOT SHAFT IN THE BEARING. 
a 
GIMBAL, (gim) 
SUBSYSTEM RECEIVING THE ANGULAR VELOCITY OF STRUCTURE CARRYING THE BEARINGS FOR THE SPINNING 
THE GIMBAL WITH RESPECT TO THE CASE AS ITS INPUT 
ROTOR OF THE GYROSCOPIC ELEMENT, ROTORS FOR THE 
AND PRODUCING AS THE OUTPUT A RETARDING TORQUE 
TORQUE GENERATOR AND SIGNAL GENERATOR, PART OF 
ACTING ON THE GIMBAL ABOUT THE OUTPUT AXIS WITH 
A MAGNITUDE PROPORTIONAL TO THE MAGNITUDE OF THE DAMPER, FLOAT SEALS AND STRUCTURE, BALANCE 
ADJUSTMENTS, STOPS, PIVOTS, ETC. 


ITS ANGULAR VELOCITY INPUT. 
RIGID ARM UNBAL ANCE COMPENSATOR, (rauc) 


SCREW ADJUSTMENTS FOR SHIFTING THE CENTER OF MASS 
OF THE GIMBAL FLOAT TO A POSITION SUBSTANTIALLY ON 
THE OUTPUT AXIS. 


* A DISCUSSION OF GENERALIZED CONVENTIONS FOR SELF-DEFINING SYMBOLS OF WHICH Al (ref}-(cmpd)}. IS AN EXAMPLE IS GIVEN BY 
DRAPER, McK AY AND LEES IN INSTRUMENT ENGINEERING (10), VOL. I. 


Fig. 2-3 Line schematic diagram fer the single-derree-of-freedenm 
integrating gyre unit. 
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The torque and signal generators are microsyn units. The rotors 
of these units are attached to the gimbal float while the stators are fixed 


to the case. 


A line schematic diagram is shown in Fig. 2-3 and a functional 
diagram in Fig. 2-4. If the gyroscopic element output is the only torque 
applied to the torque-summing member of the gyro unit, it is balanced 
by the viscous shear torque. The viscous torque is proportional to the 
angular velocity of the floated gimbal assembly relative to the case. 
This torque balance causes the output angular velocity to be proportional 
to the input angular velocity. Thus, the angle of the floated gimbal with 
respect to the case is proportional to the angular displacement of the 
case about the input axis, measured with respect to inertial space. In 
other words, the angle through which the float turns with respect to the 


case represents the time integral of the gyro input angular velocity. 


The torque generator permits introduction of command orienta- 
tion signals. The torque generator applies a torque to the float, propor- 
tional to the current supplied the torque generator. If there is no 
gyroscopic torque, the output of the torque generator causes the float 
to rotate (against the viscous shear torque of the damper) at a rate 
proportional to the current input to the torque generator. The gyro 
Output is therefore equivalent to an angular displacement about the gyro 
input axis. The output axis angular displacement is measured by the 
microsyn signal generator. The output voltage of the signal generator 


is utilized for indication or control. 


For the gyroscope to be a linear instrument, the angle through 


which the float turns with respect to the case must be kept small. For 


10 





the integrating gyro this implies closed loop operation. There are 
two types of closed (servo) loops presently in common use: 

1. a rate feedback loop in which the signal generator output 
voltage is amplified and applied to the gyro torque gener- 
ator to act as an elastic restraint, 

2.408 stabilization servo loop in which the signal generator 


output voltage is amplified and applied to a servo motor 





that supplies an angular rate about the gyro input axis. 








11 





CHAPTER 3 


INACCURACIES DUE TO SIMULTANEOUS 
VIBRATIONS ABOUT THE CASE INPUT AND SPIN AXES 


Sea General 


This chapter deals with the response of the single-degree-of- 
freedom integrating gyroscope when simultaneous angular vibrations 
are applied about the case input and spin axes, The resulting error 


torques are determined and their physical origin analyzed. 


It is assumed that when an angular vibration is present in the 
plane of the gyro case input and spin axes, components exist about both 
axes. The torque applied externally to the gyroscope is developed using 


basic gyro theory. 
3.2 Coordinate System and Symbols 
3.2.1 Coordinate System 


The coordinate systems and the relationship of the angular ve- 
locity components are shown in Fig. 3-1. The spin reference axis 
(SRA), input reference axis (IRA) and output reference axis (ORA) are 
orthogonal axes fixed to the gyro case. The gyroscopic element is 
considered to be rigidly fixed to the float so that the float input axis (IA), 


spin axis (SA) and output axis (OA) are also orthogonal. 
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A(gim) - 





SRA SA 


Warsa ~ Wacysa ©°8 Aimy + Wacyra 8" A(eim) 


Warya = Wacyra ©°S Agim) ~ Wacysa $i9 Acoim) 


Fig. 3-1. Gyroscope Coordinate System and Angular Velocity 
Components 
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3.2.2 Definition of Symbols 


The following symbols and definitions will be used in the analysis. 


A(c-r)oa = “(gim) 


IA 

IRA 

ORA = OA 
SA 


SRA 


Wgp 


WiicyIa 


WrRyIA 


Wicysa 


W(LF)SA 


Angle of float with respect to case 


Float input axis 

Case input axis 

Output axis of both the case and float 
Spin axis 


Spin reference axis (coincident with SA 


when Av = 0) 


gim) 
Angular velocity of gyro wheel about 


spin axis 


Angular velocity of case with respect to 
inertial space about case input axis 


Angular velocity of float with respect to 
inertial space about float input axis. 


Angular velocity of case with respect to 
inertial space about case spin axis 


Angular velocity of float with respect to 
inertial space about float spin axis 


Moment of inertia of wheel and float 
about float input axis 


Moment of inertia of wheel and float 
about float spin axis 
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Mapp) 


dH 

dt 
dH 
dt 











F 


M (out) 


i,j, k 


Y,6,¢ 


Moment of inertia of gyro wheel about 
float spin axis 


Total angular momentum 
Nonspin angular momentum 


Angular momentum of spinning 
wheel 


Torque applied to the gyro unit 
input axis 


Rate of change of angular momentum 
with respect to inertial space 


Rate of change of angular momentum 
in float coordinates 


Torque output about gyro output axis 


Unit vectors along orthogonal axes 
forming right-hand set 


Linear acceleration 


Linear acceleration component along 
platform axis having i as unit vector 


Linear acceleration component along 
platform axis having i as unit vector 


Linear acceleration component along 
platform axis having K as unit vector 


Linear acceleration components along 
non-rotating axes having i, j, k as unit 
vectors 


Euler angles 
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3.3 Gyro Response to Angular Vibrations 
3.3.1 Applied Torque due to Angular Vibrations 


In vector notation the torque applied\®) to the gyro is 
a. dH dH a = 
Sie as x = 
Mapp) Ff | Be se), a ae 
I F 
By Newton's Law of Action and Reaction, the gyro unit output torque is 


es aa es dH 
Mout) = ~Mappy = 2B Wary - Fa : (3-2) 


Using the relationships shown in Fig. 3-1, the following vectors may be 


written: 
War =! Worsa + J Mapa 


H = Hyg + Hgp = ig, Wagga + ly Wgp) + iy Wey) 


Werpy * #H 


EQ) -Igq) Warysa Marya * Hsp Maral (3-3) 


The relationships from Fig. 3-1 are then substituted into 


Eq. (3-3) to obtain Eq. (3-4). 


— - 2 2 
Mapp) = "LGpa - Iga) CoS” Acgimy ~ SID” A(simy) Wacyga Macy) 


2 2 
Tra ~ Iga) W'cyra — W acyga)8it Agim) ©°S A(pim) 


+ 


- Hsp (Wacyta ©°SA(gimy ~ Macysa SI" A(gimy) | (8-4) 
For a practical single-degree-of-freedom gyroscope the gimbal angle, 
A(pim) usually does not exceed one degree. (4) Therefore, the equation 
may be simplified by making use of small angle approximations and 
neglecting second order terms in Av gim)’ The torque equation is 


reduced to; 
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(ha — tga) Wacysa Macy * “Gta-Iga) Acgimy * 


2 2 
(Woicyra ~ Wacysa) ~ Asp(Wcyra - Agim) Macysa) ! 


(3-5) 
For a single-degree-of-freedom gyroscope only torques about 
the output axis are considered. The rate of change of angular momentum 
in float coordinates does not contribute an output axis torque. Equation 
(3-5) represents the applied torque of the gyro and, since it is all along 


a single axis, vector notation will not be necessary. 
a 2 2 
Mapp) =lGza ~ Iga) Wacysa Macyra * “ra ~ Iga) Acgimy“cacyra ~ “ics 


~ Agp (Wacyra ~ A(gimy Wacysa) | (3-5a) 


The terms of Eq. (3-5) represent error or undesired torques. 
A single-degree-of-freedom gyroscope, operating in a closed loop, 
would have an angular displacement from the null or zero position that 
would be in error due to these undesired torques. An examination of 


these error torques will now be made. 
3.3.2 Dynamic Torques 


The first torque component to be considered will be 


ra - Iga) Wacysa Macys (3-6) 


It is evident that the nature of the error will depend upon the variation 
of the two simultaneous angular velocities and the float inertias about 
the spin and input axes. This error is known as BeCaentee or 
dynamic unbalance. It could be eliminated by building gyro units 


having equal inertias about both axes. 
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—_—" = 2" ewe = YS 


Oe 








Next consider the second term of Eq. (3-5) 


(ra ~ Iga) Aci m) (Wacyra > Wocyga) (3-7) 
In the above expression the nature of the error will be determined by 
the difference of the inertias as well as the difference of squared angu- 
lar velocities. Variation of gimbal angle must also be considered since 


it is not a constant. 


Expressions (3-6) and (3-7) are dynamic forenoon and may be 
interpreted as representing four angular momentums and their associa- 
ted precessional angular velocities. Consider that there are four angular 


velocity components and four corresponding angular momentums. 


Momentum Components Angular Velocity Components 
Io, W 
sa “Ic)sa Wacysa 
- Agim) 41a Wc)sa - Agim) Wac)sa 
Tra Wacy1a WicyIa 
A(gim) tga Wocy1a A(gim) “(IC)IA 


The total "dynamical torque ®) when each forced precession is 


considered is: 


Mp = -Igq WacygaWacyia ~ A(gim)“acysa) 
—Tha Agim) Macysa cysa * A(gim) Wacy1a) 
+ Tha Wacyra Wacysa + A(gim) Micra) 


~ Iga Agim) Wacyra Wacyra ~ A(gim) Wacysa’ 
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The part in parentheses for each of the above terms is the forced pre- 
cessional rate while the part before the parentheses is the angular 
momentum. For example, the angular momentum component, I, AWac)s mn 
precesses about the input axis with the angular velocity, 

Wacyt bar A, gim) “(1c)s A’ When the equation is expanded and second 


order terms of A neglected, the expression for dynamic torque is 


(gim) 
obtained. 


2 2 
My = Mya-Iga) WacyraWacysa * Gra~ Iga) A(eimy(%acyra- Wacysa? 


(3-8) 
3.3.3 Gyroscopic Torque 
The applied gyroscopic torque on the gyro is 
Moe = ~HgpWacyra - Agim) Wacysa) = ~ Esp Marya 
(3-9) 


Variation of this torque will, as in previous considerations, be de- 
pendent on the angular velocities and the gimbal angle. The component, 


Hop Ay gim) Wicys A? is of special significance and is called Coning 


Error. 


3.4 Coning 
3.4.1 General Discussion of Coning 


The coning error defined above is more usually considered a 
cross coupling error. It is a result of oscillations in the plane of the 
case input and spin axes. Comparing the form of the coning expression 
from the next section with that of cross coupling reveals a close similar- 


ity. Because of this similarity, it is possible to consider cross coupling 
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as another cause of the coning error. This form of coning error may 
be reduced by restraining rotations about the output axis. A tight servo 


loop would keep the error negligible. 


If the gyro input axis has a motion such that it returns periodically 
to its original position, then the input axis will trace out a path in space. 
This motion is known as the coning motion and is present even in a theo- 


retically perfect instrument. 


The coning error is a geometric effect that allows the gyro to 
respond to an input rate produced by the coning motion. (7) Almost any 
motion of the case creates coning. Consider oscillations in the plane 
of the case input and spin axes. The component along the input axis 
will cause the gyro to precess about the output axis. When the oscillating 
motion about the spin reference axis is combined with this precession, the 
result is a coning motion of the float input axis. As another example, if 


there is oscillation of the case in the spin axis - output axis plane, coning 


would again result because the damping fluid will drag the float around. 


When single-degree-of-freedom gyros are used for stabilization, 


a coning error exists that cannot be eliminated. 
3.4.2 Platform Coning Effect °) 


Consider two coordinate systems: 


iyi, k latform coordinate system 
ty Jy p p uf 
it j PS Ky coordinate system defining the 


orientation the platform servos 
are trying to maintain 
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an do 





Fig. 3-2. Euler Angle Relations Between Platform and Desired 
Orientation Coordinates 
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By use of the Euler angle rotations,taken in the order shown in 


Fig. 3-2, the coordinate transformations, using matrix notation, are as 











follows: 
15 cos y -—sin py 0 1, 
he = sin y cos wv 0 iy 
Ky 0 0 1 Ky (3-10) 
and . 
) Wig = y Ky = alld 1 (3-11) 
i, cos 6 0 sin 9 ig 
iy = 0 1 0 Jo 
Ky — sin 6 0 cosé Ko (3-12) 
and 
w Doe cana (3-13) 
me ee al - 
iy 1 0 0 
io = 0 cos @¢ —sin @ 
| Ko 0 sin ¢ cos ¢ (3-14) 
ana = + 


It is possible to combine the three matrices and simplify them by 


making the standard small angle approximations ; 


i cosy -sinw 0|| cos9 0 sin@|{ 1 0 0 ee 
ne =|siny cosy 0 0 1 0 0 cos@ -sing I; 
TS 0 0 1j{|/-sing O cos6|| 0 sing cos ¢ EF 
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in 1 -Yy 6 Ly 
ip > V t -¢@ j 
SS —6@ @ 1 Ky (3-16) 





The angular velocity of the platform coordinate system with 





respect to the orientation coordinate system may be defined as: 


eee pit qi. ork a 
Ono pi, + 43, + 7k, (3-17) 


From Eqs. (3-10), (3-11), and (3-12) we also have 


; no = 2 * “oi Sig 
| 


® = @ es eo 
= oi, + Ojo+ Wk, (3-18) 
In order to express Eq. (3-18) in platform coordinates, substitute 


relation for ky from Eq. (3-12) and then make use of Eq. (3-14) as 


follows: 


Woo = bi, + 6 j, + PCsin 6 i, + cos 6 ky) 


= $i,+ [-P sino 6 WW cose] iy 


jg 

o 
= $i,+ [-Psin 6 cose] fl 0 0 i, 
0 cos@¢ -sin@ Jp 


0 sing cos @¢ k 


= (§ -Psin ae (6 cos ¢ + Ycos 6 sin bi, 
+ (-6 sin ¢ + ¥ cos @ cos ok, 


Assuming small angles, a comparison with Eq. (3-17) gives 
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p= ¢-0y 
a= 6 +6 

@ @ 
r= W -—d6 (3-19) 


The rates of change of the Eulerian angles are 


6 = ee oy = p+ O(r +66) = p+ré 
@ e 
= G-oy = a-¢(r+ omy = q—r¢ 
e Ed 
p= r+ = r+(q-6¥) x r+¢q (3-20) 
Integrating Eq. (3-20) to obtain the expressions for the Euler angles 
we obtain: 
t t 
g(t) = f pdt + f re dt 
O O 
t iE 
o(t)} = faqdt- f rg at 
7 O O 
f t 
Wt) = f rdt+ f a@ at (3-21) 
O O 


In the above equations, ¢, 9, and y are angles in space and when all 
are zero, the platform coordinate system coincides with the orientation 
coordinate system by definition. Having the space angles zero is no 
guarantee that the integrals of p, q and r are zero. Depending upon the 
past history of motion, p, q and r may or may not be zero. For ex- 
ample, having ¢(t) = 0 does not ensure that a @rdt is zero. A 
typical platform stabilization system can ae the integrals of p, 

q and rat zero. Therefore, the error of the platform coordinate sys- 


tem is given from Eq. (3-21) by 
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t 





A¢g(t) = f ordt 
O 
t 

Ao(t) = ~f grat 
O 
t 

) Aw) = f dqat (3-22) 

O 


If the vibratory motions are small, the integrals above become 


significant only if there is a rectification of the motions. Let 





G6 = a Sin wt 
| 0 Sn sin (wt -y) 
then 
Ag(t) = creek. r sin wt (sin wt cos y -—cos wt siny ) dt 
O 


aad 


= iA (sin? wt cos y-sin wt cos wt sin y) dt 
O 





861% wt 1 1 
= —— —— sin 2wt jcos y +—(1 + cos 2wt) sin y 
w 2 4 4 


(3223) 
An examination of Eq. (3-23) reveals one linear term. Therefore, the 


angle drifts from null with a constant drift rate. 


gd. 7 cos 7 
A 6(t) ee 


For the system just analyzed, there will exist a coning error that 
grows with time. One possible method of reducing the error would be by 
use of a coning computer. (8) The computer would solve Eqs. (3-20) and 


compensation for the coning drift could then be built into the system. 
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Seo Sculling Effect 


Sculling produces an acceleration error in accelerometers 
mounted in gyro stabilized platforms undergoing combined linear and 
angular motion. This effect will be considered because what is really 


important is the effect of this error on position indication. 


The same coordinate systems and coordinate transformation 
matrix previously used in describing platform coning apply. Define 


the total acceleration vector in both coordinate systems as 


n = nit ny Ipt By Kk, (3-24) 
n = a ie + ny Jo + 0, Ko (3-25) 
If there is no acceleration error, then 
Ree on = 0 
ny ~ "y, = 0 
ee eae = iu (3-26) 


A check on the above relations will now be made. From transformation 
(3-16) we obtain 
i = ny Up - Vip + 6 Ky] + A 


Comparison of Eqs. (3-24) and (3-26) produces the relations 
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p Xo Yo Zo 
Cae ae vi, d Bee 
La yx, Cy aa 
and by rearranging 
as "x = ae - 6 oe 
ny “ "5 = @ nze - yn 
ee Ue ey m -¢ ny. (322i) 


Equations (3-27), when compared to Eqs. (3-26), show that errors be- 
tween the accelerations actually measured and the accelerations we 


desire to measure exist. 


Assume a component of the vibration of the platform is at a 
certain frequency, w. Then consider any term of Eq. (3-27), say 
gn, - Let the linear and angular vibration vary sinusoidally at 


O 
frequency w. 


If ¢ = |¢| sin wt 
MS = Pz sin(wt - y) 
then 
é By = 19 [Pz | sin wt (sin wt - y) 
“9. 


[(1-—cos 2wt)cos y —- sin 2wt sin y | 
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which gives a steady state drift of 


#", 
cos y 





(3-28) 
As a consequence of the analysis, it is seen that small angular 
motions of the platform can result in rectified acceleration drifts whose 
magnitudes depend on the phase relationships between the linear and 
angular motions. The name "sculling" (9) is given to this effect because 
the motion causing the drift is identical to that produced in propelling a 


boat with an oar at the stern. 
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CHAPTER 4 





MAGNITUDE OF THE ANGULAR VIBRATION INACCURACIES 


4.1 General 





The errors producing forced precession of a single-degree-of- 
freedom gyroscope have been analyzed in the previous chapter. This 
chapter deals with a more practical consideration; the magnitudes of 
the various error terms. When practical,oscillations capable of being 
produced by the Angular Vibrator, built by the Instrumentation Labora- 


tory of the Massachusetts Institute of Technology, wili be assumed. 


The following are assumed as typical values for a single-degree- 


of-freedom integrating gyroscope: 


AI = I,,-Ig, = - 650 gm cm” 
2 
Hs Hap = 2x 198 gmcm" 
sec 
he | 
Cy 
A(gim) S .01745 radians 


A comparison of the various inaccuracies previously considered 


will now be made. 
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4.2 Magnitude Determination 


For convenience, the inaccuracies developed in the previous 


chapter are again listed. 


f 2 2 
Mp = Try - Iga) Wacyra Wacysa t+ “ta~ tga) Acgimycyra- Wicysa) 


(3-8) 
Mge = —Hgp Wacyta ~ A(gim) Wcysa) (3-9) 
A¢g(t) = f rx dt 
Ao(t) = - ‘ gr dt 
Ay(t) = f oq dt (3-22) 
ny, - ee = ve ~ 6 My 
yey, eh 
nM, = é Ny ~¢ By ae 


Equations (3-8) and (3-9) apply to a single gyroscope instrument while 
Eqs. (3-22) and (3-27) apply to a platform stabilized system. The mag- 
nitudes of the various terms will be determined for the case of harmonic 


angular velocities. 
4.2.1 Magnitude of the Dynamic Torque 


Let 


Wacyta = [Mya | sin (wt + 94) 


Wic)sa = sa sin (Wot + $ 9) (4-1) 
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The anisoinertia component of Eq. (3-8), when using Eq. (4-1), 
may be written 


1 
2 Al Wr Wop 





[cos [(w, — Wo)t + (9, -$5)] 





— cos [(w, + Wo)t + $,+ 99 1] (4-2) 


If 2 
W1 = Wo = 


then the above reduces to: 
Al WwW WwW [cos( ¢, — 5) - cos(2wt ) | 
: 1aA| | “SA 91 — 99) - One 
(4-3) 


Equation (4-3) consists of a constant term, and a cosine term that 
has a zero average value. Since (¢ iz $5) represents the phase dif- 
ference between the two velocities it is evident that only the in phase 
components of the motions create anerror. If the two angular veloci- 


ties are in quadrature, there will be no resultant error. 

If w, ey We , then the torque will oscillate at a beat frequency 
(w, es Wo) ; (w, = Wo) and have an average value of zero. 

The second component of Eq. (3-8) for w = W, = Wy may be 


written as: 


Al 2 2 2 
= 4 im) [ [w, | -|Wg a | + COS 2wt (|We al cos 265 


fe [Wal 2 cos 26,4) — sin 2wt (|W al 2 sin 29 9 


Ht MAN 2 sin 2¢ i)! (4-4) 
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Before the magnitude can be determined, the variation of gimbal 
angle with time must be known. For the present, an ideal integrating 


gyro will be assumed so that 


Agim) = Wacyra 


and 


A(gim) Fa f WIc)IA dt 
-|w,,l 
IA 
== COS{at+ op 
a 1 (4-5) 
Substituting (4-5) into (4-4) we get 


Si Iwral sin(wt + @,) | bal _ lWo A : 
@® 


+ cos 2wt (lWe al 2 cos 265 + Iw, al 2 cos 264) 
— sin 2ut ( [We al “ sin 2¢5 + wal 2 sin 20 »| (4-6) 


All the terms in Eq. (4-6) have zero as an average value; therefore, 
the entire expression may be neglected. The same applies for the 


case of wy yy Wo» 
4.2.2 Magnitude of the Gyroscopic Torque 
We next turn our attention to Eq. (3-9). After making substi- 


tutions for the motions we have 


M oe == Hop [Wy Al sin (w t+ 94) - A( sim) Wo A Sin (Wot + $5) | 
(4-7) 


The first term in the bracket will have a zero average value. Our 


attention will now focus on the cross coupling coning term: 
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Avgim) Hgp |Weal sin( wot +5) (4-8) 
Again assuming a perfect gyro so that relation (4-5) applies we obtain 


1 ; 
a Hop lw, IW. al cos (w,t +@ d sin (wot + $5) (4-9) 
If = b 
we can write: 


wake Hop lw, 4| Weal [sin ($5 - $4) + sin(2wt + bo + o1)] 


2W) 
(4-10) 
Thus the torque depends on the out of phase components of the veloci- 
ties. Should the two velocities be in phase there would be no torque. 
Only velocities in quadrature produce a coning error. From this 
discussion it is seen that the float drifts from null at a constant rate 


of 


1 
Std Weal sin(¢.-$,) rad/sec (Jo) 


For the condition of w, ¥ Wo, Eq. (4-9) contains no linear or 


average torque values. 


Note that when coning is a maximum, the anisoinertia effect is 
aminimum. Angular motions that are not exactly in or out of phase 
will produce both effects. However, even for this condition it is seen 
that the anisoinertia error may be neglected, since it is very small 
in comparison to the cross coupling coning error. This is due to the 
fact that the coning term involves Hop which is very much larger than 


any of the values in Eq. (4-3). 


The coning drift rate is plotted in Fig. 4-1 asa function of 


frequency for angular velocities 30, 60 and 90 degrees out of phase. 
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All of the previous work has been for an ideal gyroscope. Ac- 
tually, the true gimbal angle variation should be obtained from the gyro 
equation 


Avsimy * Wacysa A(gim) = Vaca (4-12) 


If Eq. (4-12) is used to obtain the gimbal angle motion, the only differ- 
ence from the previous result would be a nonzero average value of coning 
drift. This can be easily shown. First integrate Eq. (4-12) 


-JWacysa dt ,t I Wacysa tt 
A(gim) = © ig J’ Wacyta ee) dt 


(4-13) 
Assuming sinusoidal variations and small magnitudes for the angular 
velocities, the exponentials may be approximated by the first two terms 
of a series 


+ | W dt 


(4-14) 
Using the relations (4-1) and the gimbal angle expression (4-13), we 


can determine the nature of the coning term 





t 1 
A( sim) ft ! WA! IWa al sin (w,t + $1) cos(wot + $9) dt 
eye rae 
= —— Wal 'Woa! Pe f [[sin(w,t +6, - gt - 69] 
Wo : O 
+ sin[w t+ Wot + 61+] dt 
if Wy = Wo = 
then 


cos($ 1+ $9) - cos(2wt+, + go) 


Hl 
A(cim) =—— lA Weal sin(d, - G5) xt+ 
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It is seen that a constant drift rate appears similar to expression 
(4-11) plus an oscillatory term and a nonzero average term. Since 
the constant drift rate term dominates, the approximation for gimbal 


angle motion using Eq. (4-5) is valid. 
4.2.3 Magnitude of Platform Coning 


From Chapter 3 it was seen that platform coning effect causes 
the Euler angles to drift from null at a rate equivalent to 4 (t) 
0. Yr 


Ag(t) = % © cos y (4-15) 





The above equation was the result of integrating one of the error equa- 
tions (3-22). 
If it is assumed that 


6 
ri=er, sin (wt+ y) = W 


y=W, cos(wt + y) 
then 

vy = -W¥ow sin (wt + y) 

os -Yow 


and Eq. (4-15) is now 
-wy 6. w cos y 


AGOe= 9 a (4-16) 
Z 
The assumption will now be made that the amplitudes of the 
motions are known but that the directions and phase relationships be- 
tween the components of the motion are random. 


Let 


y K cos 8 


96 


O 
K sin £ 
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Hence 2 
ao) = 


W 





sin2B cosy 


Now assume that 6 and y are independent and uniformly distributed 


from 0 to 27, so that the mean value of each of the above trigonometric 





terms is zero, and the variance of each is 1/2. For example 


21 
Mean (cos y) = ut cos y =E(cos y) = i He cosy dy = 0 
2n 0 


2 2 a oa 
Variance (cos y) = oa (cos y) = E (cos y) = —— f cos y dy = 1/2 
27 0 


We can now write 


2 
2 
E[Ad(t) ] ia) E[ sin28 cos a 


2 (x? ‘ 1 1 
= 22 = vies 
: Coe ou 


and finally the standard deviation or RMS value is 








og = KY (4-17) 


where 
K = maximum amplitude, radians 


w = frequency, rad/sec 


Equation (4-17) is plotted in Fig. 4-2. 
4.2.4 Magnitude of Sculling Effect 


The sculling effect will be treated by the method used in the 
previous section. Consider the first of Eqs. (3-27) : 


n-n = Wn. -O6n (4-18) 
O Yo ao 
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Using the form of the sculling error found in Chapter 3 we can write 





- 1 
n-n = — n co ial 
—°o «(2 Io A ae! 2 Po", ee 
(4-19) 
If we let 
= nco 
In, | SB 
n =n sin Bp 
Ing 
| v. | =K cos a 
| 6. | =K sin a 
then 
eck 
nee m ae nK (cos a cos B cos v1, ~ Sina sin B cos Yo) 
(4-20) 


Consider that a, 6, Yp and y 9 are independent of each other and 
uniformly distributed from 0 to 27. The mean value of each trigono- 
metric function is zero and the variance of eachis 1/2. Then 


E [(n, -—n aT = ct nK)* E[ cos” a cos” B cos” et sin” a sin“ 8 cos’7 ] 
x x) 9 i 2 


1/2 
1 1 ft 1 

O(n, -n,) = —nK[ 2 (—) (—) (—) 

re 2 2 Se 
- DK standard deviation or RMS (4-20) 
4 

where 9 

n = maximum acceleration (ft/sec*) 

Ke = 


maximum amplitude (radians) 


Equation (4-20) is plotted in Fig. 4-3. 
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CHAPTER 5 


THE RESPONSE OF THE SINGLE-DEGREE-OF-FREEDOM 
INTEGRATING GYROSCOPE TO TWO ANGULAR VIBRATIONS 
SIMULTANEOUSLY APPLIED ABOUT THE SPIN REFERENCE 

AND OUTPUT AXES 


5.1 Introduction 


Gyroscopic effects caused by sinusoidal angular vibrations are studied 
in this chapter. Two vibrations of different frequencies have been applied to 
two of the body fixed axes of the gyroscope case. The particular axes to which 


these inputs are applied are the output and the spin reference axis. 


A series solution has been used to approximate the theoretical angle which 
should exist between the float and inertial space about the output axis. The 
method used was applied to a case of two equal input frequencies in ref. 10. 
Thus, this is the more general situation and a further study of the problem with 


varying inputs. 


A typical dynamic model of a gyroscope has been used as the basis for 
the equation of motion. Effects such as inertia damping, and compliance have 
been included. The resulting equation presently has no closed form solution 


and the approximate series method was used. 


The investigation was an attempt to discover if and when gyroscopic drift 


can be expected along with the associated values of magnitude. . 
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5.2 Definition of Symbols 


The following symbols and definitions have been used in the analyses: 


(units are in parentheses) 


a 


A(T_C)NS 


Ad_F)OA 


A_C)OA 


Qs FOO WD YP 


i 


constant (radian) 
constant (radian) 
constant (radian) 
damping coefficient between float and case about the 
output axis (dyne-cm-sec/radian) 
H 
“d(OA) 





Kic_F)OA 
“a(OA) 
time (seconds) 


the variable representing Ay_pyoa 


angular rotation of the case about a stationary axis 
aligned with the earth's axis which coincides with the 


spin reference axis when Ai_cyoA = 0 (radian). 
angular rotation of the float about the output axis (radian) 


angular rotation of the case about the output axis (radian) 


constant (radian) 
constant (radian) 
constant (radian) 
constant (radian) 
constant (radian) 
constant (radian) 


constant (radian) 
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H = angular momentum of the wheel (dyne-cm-sec) 


H, = constant (radian) 


lipo A = inertia of the float about the output axis (gm-cm’) 


ly = inertia of the wheel about a radial axis (gm-cm’) 
J = constant (radian) 
K = constant (radian) 

Kic_FO = elastic restraint about the output axis between float and 
case. Rotational stiffness between float and case about 
the output axis (dyne-cm/rad) 

L = constant (radian) 
= constant (radian) 
N = constant (radian) 
W, = circular frequency of vibrational input (rad/sec) 
Wo = circular frequency of vibrational input (rad/sec) 
] + 
i = ‘moa * Ww = gyroscope time constant (seconds) 
“a(OA) 


(The inertia of the wheel is added to that of the float in + 
since the wheel is assumed to be perfectly rigid with 


respect to the float about the output axis in the analyses. ) 


0.3 The Equation of Motion of the Dynamic Mode1!® 


The simplified equation of motion of the dynamic mode! shown in Fig. 5-1 
is given by : 


Kic_F)OA | Aq-moa z A-cjoa | + Ca(ga) | Agnoa 7 Sccoa | 


a 0a : Ly Aa_FIOA = 4 | Accs “(BOA | (5-1) 
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By using small angle theory, the term H | Aa cons AUL_F)O a is introduced 
by the component of angular velocity about the input axis produced by the motion 


about the space-fixed axis, Division by c d(OA) and rearrangement yields: 


K e ] + I e 
(C-F)OA F)OA*‘w]° 
ae | (roa * “qagjen * ae 2 ACL_F)OA 
d(OA) d(OA) 
Hie. | 5 q _ | 2@H0Aa 
“13 (I-C)NS “*(I-F)OA =~ A (I-C)OA * “*(I-C)OA 
d(OA) d(OA) 


or 


[ 7 Acoyns | Aa-ryoa + Aa-ryoa * 7 Aa-F0A 


= Aq_cyoa + KAq_cjoa (5-1a) 


The stiffnesses between the wheel and float about the input and output axes are 
assumed to be infinite compared to the inertia of the wheel for the frequency 
range considered. The stiffness and/or damping between the float and case 
about the input axis is also assumed infinite compared to the inertia of the 


float and wheel within the frequency range considered. 


09.4 Method of Solution of the Equations of Motion for Two Different 
Sinusoidal Vibration Inputs of Small Magnitude 


The inputs introduced to the gyroscope case may be represented as follows;2° 
A(1_C)NS = asin Wyt (522) 
A(l_C)OA = b sin Wot + ¢€ COS Wot (5a) 


In most applications of the single-degree-of-freedom integrating gyroscope, 
the constant h is designed equal to unity. For purposes of simplification, this 


constant will be included in the vibration inputs as a factor of the constant a. 
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To be strictly correct, earth rate should be included in the definition of the 
vibration input about the North-South axis. This rate is, however, small 


compared to the input angular velocity. 


Substitution of equations (5-2) and (5-3) into equation (5-1a) yields a 


linear differential equation with a variable coefficient of the following form: 


ee 


| k — aw, cos ot | X + x + Ts 


- bwo COS Wot — CWo sin Wot + kb sin wot + ke cos Wot (5-4) 
where AU_F)OA = X 


There is no known closed form solution for equation (5-4). However, if the 
constants are small compared to unity, a suitable series solution may be ob- 


tained. * 


5.4.1 Negligible Elastic Restraint Between the Float and Case 
About the Output Axis 


When the elastic restraint is very small compared to the damping con- 


stant, i.e., 


k>0 
equation (5-4) reduces to: 
— aw, COs wt) x+ X+ 7X = bwo cos Wot — CWo sin Wot (5-5) 
Assume a solution of the form: 
X = X; + Xp (5-6) 
where xy is defined by: 
Xy 2h TK = bw» COS Wot — CWo sin Wot (52a) 


Substitute (5-6) and (5-7) into (5-5): 


(-aw, cos ws t)(xy Et Xo) A (x, + Xp) a 7(X, Ee Xp) = bwo cos Wot - CWo sin 2 | 
5-8 


*The series solution technique was suggested by Mr. Daniel Goldenberg 
of the Mathematics Group, MIT Instrumentation Laboratory. 


ALC 





or 


(-aw, cos Wt) xX» + X> + Tos = (aw, cos wt) Xy (5-8a) 
Now let x, = Xe + Xy (5-9) 


where X9 is defined by: 


se + ie - (aw, cos wt) xy (5-10) 
Substitute (5-9) and (5-10) into (5-8a): 
(~aw, cos Wt) (xX, - X4) + (X, te x,) a T (Xz + X,) = (aw, cos wt) Xy (5-11) 


or 


(-aw, cos wt) X4 + a + TK, = (aw, cos wt) Xq (5-11a) 
Continuation of this process results in a series solution of the following form: 
X= Xp + Xgt Reto... eee Xo 1-1 + Xam (5-12) 


where Xom is the even term remaining after the above process is repeated 


m times. Xm satisfies the equation: 
(—aw, cos wt) yey Gree US yen a (aw, cos wt) Xo] (Seas) 


A Brief Summary for k> 0 is: 


Problem: 
(—aw, cos w,t) x + eee dW, COS Wot — CWo Sin Wot (525) 
soe eee Xe Xe te we OX + X (5-12) 
1 3 ° 2m-1 2m 

where: 
fe + 7a = bw, COS Wot — Cwo Sin Wot (5-7) 
a + axe = (aw, cos wt) Xy (5-10) 
X. + ik = (aw, cos wt) Xe (A-6) 


X(one1) + 7 *(2N+1) = (aw, cos wt) X(2N+1) (A-6a) 
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Each successive term in the series of equation (5-12) can be shown to be 
smaller than the preceding term by a factor of approximately a. (Refer to 


3 radians, 


ref. 10 for further discussion.) Since a is of the order of 2 x 10° 
this series is seen to converge rapidly. Hence all terms beyond X. can be 


reasonably neglected. 


The following evaluation of the series and associated coefficients clearly 
shows the appearance of the factor a which permits terms beyond a certain 


point to be neglected. (See Appendix A for solution details. ) 


Ad_F)OA = x k = elastic restraint = 0 

X = X, + X3 + Xp + other terms which are neglected (5-12) 

X, = A sin Wot + Bcos Wot (A-7) 

X_ = C sin Wat + D sin Ww gt + E cos Wot + F cos wyt (A-17) 
We = Wy + Wo Wa = Wy - Wo 


Xx, = G sin (2w,t + Wot) + Hcos (2w,t + Wot) + Jsin (2w,t - wo) 


+ Kcos (2w,t — wot) a L sin (—wot) + M cos (—wot) 
+ Nsin oD + P cos (wot) (A-37) 
b+cTw C-bT Wo ane 
See poe - 
eee ee) 1 a, 
Ler Wo + T Wo 
aw,(B —7T Aw.) aw,.(B+ 7TAw,) 
eee Ss (A.-228) p- __!._* Ea 
W(l + T° Ws wall + 7" wy 
—aw,(A+ TBw,) aw.(A —7Bw,) 
pee (A228) F = on (A-30) 
2W9(1 i 708) 2w4(1 a Ww) 
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aw, [E- 7C(2w, + wo) | 


~aw, [C + TE(2w, + Wo) | 


en = 
2(2w, + Wo) [1 + 7*(20, : wo)” J 2(2w, + Wo) [1 + 7*(20, + wo)” | 
(A-4'7) (A-48) 
= AW, [F - TD(2w, ~ Wo) ] — —aw, [D + TF(2w, = Wo) | 
2s] = Wp) [1+ 77(2w, - 0)" ] 2(2w, - 5) [1 + 72(2w, - w9)? ] 
(A -50) (A-49) 


aw, [-E + TCw, | 


Zee 
2Wo [1+ TW | 


- C E 
(A-51) M = aoe ee (A-52) 


2 2 
2Wo [1+ 7° we | 


L= 


aw, [F + TDw, | aw,(D —T Fw, | 
eee (A-53) D . .a cae (A-54) 


Zee 2,2 
2Wo [1+ 7” wo | 2Wo [1+ 7° wo | 


.0009 sec = .0 X 107 3sec 

3 Typical values as used by 
radians the MIT Instrumentation 

Laboratory Angular Shaker 


T 


Amax 2x 10° 
qib2+c2)_ = 4x 10-2 radians 


With negligible elastic restraint between the float and case, the series 
solution consists of sinusoids. Such sinusoids will exhibit a time-average 


value of zero except for special cases which are of interest. 
The first case has been studied in ref. 10 and concerns equal frequencies, 
Wy = Wo 
As in further explained in Appendix A, equal input frequencies will cause the 


terms of D sin wat and F cos w,t to have indeterminate values. By applying 
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L' Hospital's Rule, a constant drift term which is time dependent is found 
to exist. This term is: 


Drift = | t where w = W1 = Wo 
2 





This existing drift is identical to that found in ref. 10. 


The second case of interest concerns the almost natural region for sus- 
picion, namely harmonics. Terms concerning Xp of the series give indetermi - 
nate values when 204 = Wo. By L' Hospital's Rule, a drift term appears which 


is very similar to the term arising from equal frequency input. This term is: 





Faw, 
Driit = t See Appendix A for proof and 
2 necessary mathematics 


204 = Wo 
If further series terms such as Xr, and Xg were evaluated, drift terms 
associated with other harmonic multiples are found. Thus, the conclusion is 
that equal frequencies and harmonics are capable of causing drift for the case 
of negligible elastic restraint. 
The physical case wil! probably display many superimposed input fre- 
quencies. This is likely to have some harmonics thereby creating drift. Such 


drift is substantiated by this theoretical solution but it is suspected that fre- 


quencies approaching the harmonics may well create actual physical drift. 


5. 4. 2 Finite Elastic Restraint Between the Float and Case 
About the Output Axis 


For the case when k 7 0, it is desired to find the solution of the 


following equation: 


(k-aw, cos wt) x 4X4 7K = bwWo cos Wot — CW sin Wot + kb sin wot +kccos Wot 


(5-4) 


00 





The mathematical solution with details is in Appendix A. The resulting series 


with values of coefficients in terms of the known constants is as follows: 





X= X; + Xg +X, + other terms which are neglected (5292) 
| 
, x, = Asinwot + Bcos Wot (A-67) 


XQ = D sin w,t sin Wot + E sin w,t cos Wot 





+ F cos w,t sin Wot + G cos wt cos Wot (ASH) 


Xe = H, sin 2w,t sin Wot + I sin 2w,t cos Wot 
+ J cos 2w,t sin Wot + L cos 2u,t cos Wot 


+ M sin wot + N cos Wot (A-89) 


Expressions for the coefficients in terms of known constants may be found in 


Appendix A. The coefficients have been left in determinant form. 


By the use of suitable trigonometry relations the following forms are 


obtained for XQ: 


Xg = = [cos(w, —wo)t] - = [cos(w, + wo)t] 
i ~ [sin(w, + Wo)t] + a [sin(w, _ Wo)t] 
+ e [sin(w,+w)t] - = [sin(w, — Wo)t] 
G G 
see t — = t A-15 
1 ; [cos(w,+Wo)t] + 3 [cos(w, — w.)t] ( a) 


When W1 = Wo, the case of ref. 10 occurs. An indeterminate form does not 
appear because k is finite. Thus time dependent drifting is not a factor with 
non-zero elastic restraint. However, this special case of equal frequencies 
gives a value of unity to terms containing cos(w, ~ Wo)t and the result is a 
steady non-zero value for the angle between the float and inertial space about 


the output axis. This is verified by ref. 10. 
o1 











The other case of interest with finite elastic restraint again concerns 
harmonics. This can easily be seen by converting Xp terms by means of 


trigonometry relations to the following: 


H H 
Xp = = [cos(2w, -wo)t] - — [cos(2w, + w»)t] 
+ es [sin(2w,+o)t] + - [sin(2w, - w,)t] 
+ = [sin(2w,+w,)t] - + [sin(2w, - w»)t] 
L L ~ 
+ —— [cos(2w,+w.5)t] + — [cos(2w, -—w.,)t 
: 1+ %o)t] ; [ By ie ] 
+ M sin Wot + Ncos Wot (A-89a) 


The drift resulting from the harmonic when 2W4 = Wo is clearly seen as 
a possibility. However, because ofa finite value for elastic restraint, indeter- 
minate forms do not appear. Nevertheless, this particular condition causes an 
angle to take place about the output axis between the float and inertial space. 


This occurs when cos(2w, - Wo)t is unity in value. 


This case gives a term showing that a constant angular displacement 


occurs for the harmonic condition instead of a time dependent drift term. 


Such a condition is almost a natural situation for a gyroscope. As the 
finite restraint is increased, the displacement between the float and the refer- 
ence in inertial space should become smaller. The extreme limits are for 
zero elastic restraint and infinite elastic restraint. The former case gives a 
time dependent drift while k =o creates a situation permitting no angle, what- 


soever, to exist between the float and inertial reference. 


The most interest undoubtedly lies with the case of finite, non-zero, 


elastic restraint. This is the existing case of practical importance. The 
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concluding results are that 1) a time-independent drift term is possible if 
input frequencies contain harmonics or equal values, 2) float angular dis- 
placement with reference to inertial space depends on the degree of elastic 
restraint in the gyroscope, and 3) high harmonic terms may possibly be 
neglected because of their small value which becomes insignificant when 


compared to manufacturing tolerances. 


It is pointed out that a drift may appear to be present at any one instant 
of time if a time-averaged measurement is not obtained even if harmonics are 
not present. This could be of noticeable magnitude dependent upon the input 


frequencies and their own absolute magnitudes. 


It is strongly recommended that this theory be applied to the MIT Instru- 
mentation Laboratory angular shaker in an attempt to verify the results. Har- 
monic inputs appear feasibile for the shaker by minor adaptation after sufficient 


data has been collected concerning the case of two inputs of equal frequency. 


Terms in the approximate series solution beyond Xp have only theoretical 
interest at present. Until gyroscopic tolerances exceed their present values, 
common sense tells us that the selected series is as good for analysis as would 
be the unknown closed form solution for the equation of motion. 

5.9 Drift of a Typical Gyroscope for the Case of Negligible Elastic 

Restraint About the Output Axis with Two Equal Input Frequencies 

The ensuing results are based on the series solution for the equation of 
motion when the elastic restraint between the case and float about the output 
axis is negligible. A further stipulation is that w, = W» in order that this in- 
formation may apply directly to the MIT Instrumentation Laboratory angular 


shaker. 
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The angular shaker has a frequency range from twenty (20) cycles per 
second to one hundred (100) cycles per second. The constant a" has a max- 
imum value of .002 radians and the constants '"'b'' and ''c"’ are limited such 
that b2 * oe has a maximum of .004 radians. A typical value of 7, the 
time constant of the damper-mass system, is .0009 seconds. 

By definition: 
Ad_cyns = a sin w,t (5-2) 


AL_cyoa = bsin Wot + ¢ CoS Wot (5-3) 


For this particular case, w = W1 = Wo as the angular vibrator is capable of 


applying only one frequency at a time to the two input axes. 


ienca—n0) A(r_c)NS and Ai_cyo A are in phase and they will differ only 
by the magnitudes of "a" and "b". (a = bis a special case.) In general, "'b” 
and ''c”’ will have finite values thus causing A(I_c)Ns and A_c)o in to differ 
by a phase angle as well as magnitude. (Here, again, a=.{b2+c* isa 


special case. ) 


7 1 2 2 
b sin wt + c cos wt ofa + cot 1b/c } b +c (5-14) 


cot b/e = @, a phase angle 


If W1 = Wo, the following is drift rate for k 0: 


1-— Tw 
Drift Rate = 26%} __¢ (5-15) 


{ee 


A plot is made for drift rate /ac for various ratios of b/c. 


es TW Pere 

Drift Rate _ ee 2 |S ee (5-16) 
— = 12 2 2 2 2 2 
ac 1+7 @ 1 42% 2G) 
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Positive and negative values of b/c were selected to correspond to phase angles 
from approximately + 5° to + 84°. Ratios for b/c selected are indicated on the 
graph. 

The curves of Fig. 5-2 show that for a frequency range of 20 to 100 cycles 
per second, the drift rate will undergo a change in direction for b/c = + 2.0 
around 90 cycles per second. The frequency required for this reversal is fur- 


ther decreased for increasing values of b/c. 


Although the curves represent specific values for constants, it is believed 
to be typical of gyroscopes which may be analyzed by the MIT angular shaker. 
It is suggested that these curves be compared to actual data from the vibrator 


test machine at the earliest opportunity. 
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CHAPTER 6 


CONCLUSIONS AND RECOMMENDATIONS 


6.1 Conclusions 


1. When a floated single-degree-of-freedom integrating gyro- 
scope is subjected to vibrations undesired torques or drift 
rates may exist. A summary of such unwanted errors, 
arising when the vibrations are angular in nature, is pre- 
sented in table form. The errors are classified by name, 


type and cause. 


Name Type Error Cause 


Anisoinertia effect Dynamic Simultaneous angular 
motions about the case 
input and spin axes 
along with unequal mo- 
ments of inertia about 
these axes. 


Cross coupling coning Gyroscopic Float precession about 
its output axis which 
causes a component of 
the case spin axis angular 
velocity to be coupled to 
the float input axis angu- 
lar velocity, 
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Platform coning or 
gimbal walk 


Sculling effect 


Geometric 


Geometric 


Simultaneous angular 
velocities about plat- 
form axes such that 
the float input axis 
undergoes a periodic 
conical motion 


Simultaneous angular 
motion about any plat- 
form axis and linear 
acceleration along a 
perpendicular axis 
producing a rectified 
acceleration along an 
axis perpendicular to 
both axes. 


2. A second table is presented to illustrate which effects are 


present when a single-degree-of-freedom integrating gyro- 


scope is subjected to constant angular velocities and periodic 


vibrations which may be in or out of phase. 


Constant Angular 


Velocities 
Anisoinertia Yes 
Cross coupling Yes 


coning 


Platform coning Yes 
or gimbal walk 


Sculling* Yes 


Periodic In Phase 
Angular Velocities 


Periodic Out of Phase 
Angular Velocities 


No 
Yes 


Yes 


No 


*For sculling the phase refers to the phase relation between 


the linear acceleration and the Euler angle motion (space angle motion). 
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Coning is a cause of drift for a single-degree-of-freedom 


integrating gyroscope. 


. Cross coupling coning is inversely proportional to frequency 
and directly proportional to the product of the amplitudes of 


the angular motions. 


. Cross coupling coning may be minimized by use of a tight 


servo loop which keeps the gimbal angle close to zero. 


. Platform coning is directly proportional to frequency and di- 
rectly proportional to the square of the amplitude of angular 


motion. 


. Sculling generated drift is directly proportional to the ampli- 
tude of the angular motion and linear acceleration. 


. For zero elastic restraint between the float and case about 
the output axis: (Based on a series solution of a dynamic model) 
(a) Constant drift rate appears when two equal frequencies 
are simultaneously applied to the spin and output axes 

(e.g., wy = Wo ) 

(b) Constant drift rate appears when the frequencies applied 
simultaneously to the spin and output axes have a har- 
monic relation (e.g., 2w, = Wo) 

(c) The total drift at any instant of time consists of several 
superimposed sinusoids. Except for the special cases 
concerning inputs of equal or harmonic frequencies re- 
ferred to in (a) and (b), the time averaged value of 
angular displacement of the float with respect to inertial 


space about the output axis is zero. 
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9. For finite elastic restraint between the float and case 
about the output axis: (Based on a series solution of a 


dynamic model) 





(a) Drift terms which are time independent appear 
when two equal frequencies are simultaneously 
applied to the spin and output axes (e.g., W, = Wo). 
These drift terms have a constant magnitude de- 

| pendent upon a, b, c, 7, Wy; and Wo. but not 

dependent upon time. 

(b) Drift terms which are time independent appear 
when the frequencies applied simultaneously to 
the spin and output axes have a harmonic relation 
(e.g., 2w, = Wo). These drift terms have a con- 
Stant magnitude dependent upon a, b, c, 7, W4> 
and Wo but not dependent upon time. 

(c) The total drift at any instant of time consists of 
several superimposed sinusoids. Except for the 
special cases concerning inputs of equal or har- 
monic frequencies referred to in (a) and (b), the 
time averaged value of angular displacement of the 
float with respect to inertial space about the output 


axis is zero. 
6.2 Recommendations 


1. Verify theoretical results by use of the MIT precision angular 
vibrator. This should be done for various values of elastic 


restraint. 
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. The present vibrator is capable of applying two simultaneous 
inputs of equal frequency to two orthogonal gyroscope axes. 

It is suggested that the machine be investigated to determine 
if input frequencies with harmonic relationships can be simul- 


taneously applied. 


. Linearity should permit results obtained from a specific fre- 
quency applied to one axis to be superimposed upon the results 
obtained from a different frequency applied to another gyroscope 
axis. With this procedure, it may be possible to obtain physical 
data concerning simultaneous inputs of any desired frequencies 
within the limit of the vibrator. Information for assorted com - 
binations of input frequencies to two orthogonal gyroscope axes 
could thus be obtained from an angular vibrator capable of 


applying only one frequency at a time. 


. Great care should be taken to obtain valid data in the low fre- 


quency range. This is the suspected source of many uncertainties 
and may be caused by the relatively obscure bending modes of 


large missiles and high speed aircraft. 
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APPENDIX A 


SERIES SOLUTION FOR THE EQUATION OF MOTION OF A 
SINGLE -DEGREE-OF-FREEDOM INTEGRATING GYROSCOPE 


A.1 Series Solution for Negligible Elastic Restraint Between 
the Float and Case About the Output Axis (k->0) 





The general equation is: 
(k -aw, cos wt) X4+X4 7K = bwWo cos Wot — CWo sin oot kb sin Wot 
+ ke cos Wot (A-1) 
If k>0, we have: 
(-aw, cos wt) X4+X4+ TX = bwo cos Wot - CWo sin Wot (A-2) 
In accordance with the series established in Section 5.5, the interest lies with 


the following terms: 


X= Xy + Xg + Xp tee ee + XQ 4 + Xam (A-3) 
where: 
X4 + TK, = bwWo cos Wot — CWo sin Wot (A-4) 
Xq + TX. = (aw, cos wt) Xy (A-5) 
Xp + TX, = (aw, cos wt) Xe (A-6) 
(ON 1) + TX ony 1) ~ (aw, cos wt) X(2N+1) (A-6a) 
Let 
ee = A sin Wot + Bcos Wot (A-7) 
1 = A Wo cos wot - B wo sin wot (A-7a) 
25 = -~A we sin wot - B we COS Wot (A-7b) 
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Substitute (A-7), (A-7a), and (A-7b) into (A-4) 


[Awe cos wot- Bwo sin wot] + 7[-A ws Sin Wot - Bws COS Wot] 
= bwWo cos Wot — CWo sin Wot 
Equate coefficients of like terms and solve for A and B. 


AW, -TB ws = bwo (from cosine terms) 


—Bwo - TA ws = =CWo (from sine terms) 


A= b+ TB Wo 


Be C—- TAW, 


A = b+ Two[c-TAws] 


A + Two A =e CT Wo 





B = c-Tw,[b+ TBw,] 


2 ae, _ 
B+T Wo B - C-—bT Wo 





Only those terms affecting the steady state solution are considered. 


An equation for x, is obtained by substituting (A-7) into (A-5). 


X, + TX, = (aw, cos w t(A sin Wot + Bcos Wot) 


"s + TX, = Aaw, cos wt sin Wot + Baw, cos w,t cos Wot 
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(A-8) 


(A-9) 


(A-10) 


(A-9a) 
(A-10a) 


(A-11) 


(A-12) 


(A-13) 


(A-14) 


(A-15) 


(A-15a) 





From basic trigonometry relations, the following are obtained: 


cos w,tsinws.t = [sin(w,t + Wot) _ sin(w,t = Wot) ] 


as 
2 
4 





cos w,t cos wot = [cos(w,t + wot) + cos(w,t - wot) | 
Let wet = (wt + Wot) and w,t = (wt _ Wot) 


In view of these relations: 








ee Aaw Baw 
Xa t+ TX. = [sin wet — sin w,t] + [cos wot + cos wyt] (A-16) 
Let 
X_ = C sin Wat + D sin wat + E cos Wat + F cos wat (A-17) 
Ss = Cwe cos Wet + Dw, cos Wat - Ewe sin Wet — Fw, sin wat (A-17a) 


- Cws sin wat -Dw4 sin Wat - Ews cos Wel - Fw4 cos wat (A-17b) 


Substitute (A-17), (A-17a), and (A-17b) into (A-16) . 


[Cw, cos wat + Dw, cos wyt - Ew, sin wet — Fwy sin wt] 


+ [-7Cws sin wot -TDw4 sin wat -TEws COS Wot -T Fw4 cos wat | 


Aaw Aaw Baw Baw 


sin Wet - sin wat + 5 cos Wat + 1 cos Wat (A-18) 














Equate coefficients of like terms and solve for C, D, E, and F. 











9 Baw, 

COS wat: Cw. -TEW, = ; (A-19) 
9 Baw, 

cos wat: Dw, -TFwy = . (A -20) 
9 Aaw, 

sin Wet : -Ew. -TCws = ; (A-21) 
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Aaw, 















sin Wat: -Fw, -1Dw4 = = ; (A-22) 
1/2 Baw, + TEWS 
| _ (A-19a) 
a: 
) 1/2 Aaw, - 7Cw 
——————— (A-21a) 
“3 
| Substitute (A-21a) into (A-19) 
9 | -1/2 Aaw, -7Cws Baw, 
w 2 
3 
Cw. + 1/2 TAaw Ws +7Cw3 = 1/2 Baw, 
C(w. + 7? w) = 1/2 Baw, —1/2 rAaw,w 
3 3 1 3 
Baw, — TAaw, w 
2(we + 77 w3) 
Substitute (A-19a) into (A-21) 
9 1/2 Baw, + TEWS Aaw, 
EO TW ee ee (A-25) 
3 3 9 
<3 
Sho, —1/27Baw, wo _-EW? = 1/2 Aaw 
3 3 3 1 
E(we f r* ws) = = f2 Aaw, Sei. T Baw, Ws 
- Aaw, - TBaw,w 





2(w4 +77 w) 
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2 
1/2 Baw, + TF w, 





4 
1/2 Aaw, - TD we 
se 1 4 
4 
Substitute (A-22a) into (A-20) 
9 1/2 Aaw, - TD w4 Baw, 


Dw, - TW POU as =! tate 
4 4 W4 9 








Dw, = 2 TAQW Wy + r*Dws = 1/2 Baw, 
Diw, + 77w°) = 1/2 Baw, + 1/2 rAaw,w 
4 4 1 14 


Baw, + TAAaw, Wy 





2(wW4 + 7 we) 


Substitute (A-20a) into (A-22) 


9 1/2Baw, aL T Fw, 


-Fw, - Tw, | ——___—__—__| = = W/2 Aaw, 
aA 
Fw, -1/2 rBaw,w, - T*Fw? = — 1/2 Aaw 
ed Boer tend 4 1 
F(w, + 7% w) = 1/2 Aaw, -— 1/2 TBaw,w 
4 qe 1 1°4 


Aaw, a TBaw, wy 





2(w4 + 7 4) 


An equation for x, is obtained by substituting (A-17) into the following: 


Xp + T%p = (aw, cos 4) Xe 


(A-20a) 


(A-22a) 


(A-27) 


(A-28) 


(A -29) 


(A-30) 


(A-31) 


oe + TX. = aw, cos wt [C sin wot + D sin wt + E cos wet + F cos wat] (A-32) 
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oe e 
TXe+X— = Caw, cos wyt sin Wet + Daw, cos wt sin wat + Eaw, cos wt cos Wet 


+ Faw, cos w,t cos wat 


Trigonometry Relations 


cos w,t sin Wet = ~ [sin(w,t + Wat) - sin(w,t = Wet) | 
cos wt sin wat = = [sin(w,t + wt) - sin(w,t _ wat) | 
cos wt cos wet = - [cos(w t+ wot) + cos(w,t — wet) ] 
cos w,t cos wyt = = [cos(w,t + wyt) + cos(w,t — w,t) ] 
Wet = (wit + Wot) 
wat = (wt -— wot) 


1 , ‘ 
TX, + Xp = oI Caw, [sin(2w,t + Wot) _ sin(—wot) | 


p= 


a Daw, [sin(2w,t —wot) — sin(wot) | 


bh 


+ a Eaw, [cos(2w,t + Wot) a cos(—woft) ] 


1 
os Faw, [cos(2w,t —wot) + cos(wot) | 


(A-32a) 


(A-33a) 
(A-33b) 
(A-33c) 


(A-33d) 


(A-34) 


(A-35) 


(A- 36) 


(A-37) 


Let 
X, = G sin(2w,t + Wot) + H, cos(2w,t + Wot) | sin(2w,t - Wot) + K cos(2w,t — Wot) 
+ L sin(—w,t) + M cos(—Wot) +N sin(w, t) + P coS(Wot) 
a = G(2w, + Wo) cos(2w,t 7 Wot) = H,(2w, 4 Wo) sin(2w,t + Wot) 


+ 
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J(2w, _ Wo) cos(2w,t = Wot) “i K(2w, = Wo) sin(2w,t = Wot) 


Lwo cos(—wot) + Mw» sin(—Wot) + Nw, cos(wot) ze Pwo sin(wot) 


(A-37a) 





Xp = - G20, + we)” sin(2w,t + Wot) = H,(2w, + wy)? cos(2w,t a Wot) 
28 
-JS(2w, _ Wo) sin(2w,t - Wot) = K(2w, _ wo)” cos(2w,t = Wot) 


= Lws sin(-wot) — an ws cos(— wot) = Nws sin(w ot) —P ws cos(wot) 
(A-37b) 
Substitute (A-37), (A-37a) and (A-37b) into (A-36). 


-TG(2w, + wo)? sin(2w ,t + Wot) —T H (20, + wo)” cos(2w,t +. Wot) 

-TI(2w, _ Wy)” sin(2w,t _ Wot) = TK(2w, = toy)" cos(2w,t = Wot) 

Tels we sin(—wot) - T™M ws cos(—Wwot) = Ti Nws sin(wot) = Pws COS(Wot) 

+ G(2w, + Wo) cos(2wyt + Wot) _ H,(2w, + Wo) sin(2w,t + Wot) 

- J(2w, a Wo) cos (2w,t = Wot) a K(2w, = Wo) sin (2 ,t- Wot) -Lw» cos(—Wwot) 
+ Mw, sin (—wot) + Nwo cos(wot) - Pwo sin(wot) = 2 Caw, sin(2w, + wot) 
= 1/2 Caw, sin(—wot) 4 1/2 Daw, sin(2w,t —wot) mole Daw, sin(wot) 

1/2 Eaw, cos(2w,t + Wot) + 72 Eaw, cos(—Wwot) lye Faw, cos(2w,t — Wot) 


at / 2 Faw, cos(wot) (A-38) 


Now equate like parts and solve for coefficients. 


sin(2w,t+ wet): -TG(2w, + W»)*- Hy (20, +0») = = Caw, (A-39) 
cos (2w,t “ Wot): -7TH,(2w, + Wy)” + G(2w, + Wo) = = Eaw, (A-40) 
sin (2w,t = Wot): -TI(2w, — Wo)” ee K(2w, = Wo) = - Daw, (A-41) 
cos (2w,t = Wot): -TK(2w, = Wo)? i J(2w, _ Wo) ~ * Faw, (A-42) 
sin(—wo1): -rLws + Mw, = - 1/2Caw, (A-43) 
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: cos (- wot) : —-™M ws - Lw, = + Baw, (A-44) 


eo 2 
sin (wot) : -TNw, — Pw, = -= Daw, (A45) 
cos (wot) : —-TP we + Nw, = i Faw (A-46) 
a 2 2 | 9 1 


- 1/2Caw, - TG(2w, + wo)? 
(2w, A Wo) 


Substitute (A-39a) into (A-40) and solve for G: 


= Ve Caw, —T G(2w, + wo)? 


2 
ane (2w, ae Wo) + G(2w, + Wo) = = Baw, 


(20, + Wo) aed 


1 sa ee 3 aot 
ce TCaw, (2w, + Wo) + T. G(2w, + Wo) + G(2w, + Wo) Se Eaw, 


G (20 + Wo) - 7 (2w, “ 0 | = = Eaw, = TCaw, (20, + Wo) 
Eaw, = TCaw, (2w, + Wo) 


G = (A-47) 
2[(2w, + Wo) + 7 (2w, + We)? | 





or 


aw, [E— TC(2w, + Wo) | (aan 


2[(2w, + Wa) + 77 (2w, + Wy)? ] 





2 
1/2 Eaw, + T H,(2w, + Wo) 
pee (A -40a) 
(2w, + Wo) 


Substitute (A-40a) into (A-39) and solve for H, : 


2 
1/2 Eaw, + 7 H,(2w, + wo) 
-T (20, + Wy)° ooo! he eee = H,(2w, + Wo) -1 Caw, 
(2w, + Wo) 2 
(A-39) 
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myi 3 > A 

pee 1s * 2) — TH, (2w, + wo) — H, (2H, + wo) = es Caw, 
H, [(2w + Wo) + 7" (2w oe 1 Ca 1 + Ea (2 ) 
1 fr 7 5 et 5 SS = mee, + 8 


~Caw, -7 Baw, (20, + @,) 
1 (A-48) 


2[(2w, + Wo) + r(2w, + wo)? 
or 


aw, [-C -T E(2w, + Wo) | 
ES a (A -48) 


2[(2w1 + Wo) + 7 (204 + wo)? | 





1/2 Faw, + TK(2w, -w)" 
J = = ae (A-42a) 
(20, os Wo) 


Substitute (A-42a) into (A-41) and solve for K: 


2 
\2 1/2 Faw, + TK(2w, x Wo) 


-1 (21 - Ws -K(2w, - 9) = = Daw, (A-41) 


(2, — Wo) 
-~+ 7 Faw, (2w — Wo) — 7 °K (2w -w,)° - K(2w,-@) ==> Daw 
[ imei “2 Ieee pe "9 l 


1 1 
K[ (20, = Wo) + 7 *(2w = wo)? | = — = Daw, - a 7 Faw (20, - Wo) 









AW, [-D-T F(2w, = Wo) ] 


K = (A-49) 
2[(2W4 - Wo) rt 77 (2w, as Wy)? ] 
2 
Sie2 Daw, - TI (20, = Wo) 
Ke OEE EEE (A-41a) 
(20, _ Wo) 


710 








Substitute (A-41a) into (A-42) and solve for J: 


2 
=1/2Da0. —-7Tal2w, —Wa) 
-7 (2, - #9)” ee eee 


1 2 3 =! 
a i Daw (2w, -_ Wo) at J(2w, - Wo) + J(2w, - Wo) A Faw, 


J[ (2w, - Wo) + 7 (20 - wo)? | = = Faw, - = TDaw, (20, _ Wo) 









aw, [F- TD(2w, - w») | 
2 3 
= 1/2 Cag me Ww? 
1 Z 
2 
Substitute (A-43a) into (A-44) and solve for L: 


o| - 1/2Caw, ern ws 


SA 
=T Wo | ——____—_—_ -Lw, =o Eaw, 
Wo 2 
1 2 3 = ol 
oo Caw, Wo ih Lw = Lwo = a Baw, 
L(we - 7? we) = ne Eaw, a T Caw, Wo 





2 
= 4/2 Eaw, -TMw5 


3 
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1 
* J(2w, - Wo) = 3 Faw, 


(A-42) 


(A-50) 


(A-43a) 


(A-44) 


(A-51) 


(A-44a) 





Substitute (A-44a) into (A-43) and solve for M: 


9 -1/2Eaw, -TMws 1 
Se Ne a eee 
Z 2 1 

Wo 2 


1 2 3 _ 1 


M(w» + r ws) = =i Caw, -< T Baw, Wo 





2 
; 1/2Daw, = TNw, 


ee) 
Substitute (A-45a) into (A-46) and solve for N: 


9 1/2Daw, -TNws 1 
-T Wy | | + Nwg = = Faw, 


N(wo a. 7 ws) = “ Faw, a T Daw, Wo 





AW, (F + T Dw) 


2(Wo + 7 we) 


1/2 Faw, +7P ws 


5) 
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(A-43) 


(A-52) 


(A-45a) 


(A-46) 


(A-53) 


(A-46a) 








Substitute (A-46a) into (A-45) and solve for P: 


9 1/2 Faw, + 7 Pw 


io ee = ee, (A-45) 
Wo 2 
1 2 3 1 
es = Fi Faw Wo —T Pw _ Pwo = ae Daw, 
P(Wo + rw) = = Daw, -— t Faw, Wo 


AW (D-7T Fwo) 
ae (A-54) 


2(wo + rw) 





X= X, + Xgt+ Xp +... - Xp 4 + Xp (A-3) 
x, = A sin Wot + Bcos Wot (A-7) 
Xg = C sin Wot + D sin wat + Ecos Wet + F cos wyt (A-17) 


Xp = G sin (2w,t + Wot) + H, cos (2w,t + Wot) + J sin(2w,t - Wot) + K cos (2w,t — Wo! 


+ L sin (—wot) + Mcos (—wot) + N sin (wot) + Pcos (wot) (A-37) 


Terms beyond xX, may be neglected if the order of magnitude of a is much 
smaller than unity (maximum a on the MIT angular shaker is 2 x joe 
radian). If xy is defined to be of the order of magnitude of one, Xg is of the 
order of magnitude of a, and Xp is of the order of magnitude of a’. Hence 
each successive term in the series of x is smaller than its predecessor in 


accordance with the value of a. 
For further discussion on convergence and a process for iteration, 


ref. 10 is recommended. 


73 











If only the terms associated with X1> Xg, and X, are considered, it is 
readily apparent that the time average of drift for the sinusoids is zero. Close 
scrutiny, however, shows drift possibilities as Wy and Wo approach one another 


in value and also when harmonics are present. 


The specific case where W1 = We has been investigated in ref. 10. The 
general case which is the subject of this chapter degenerates to the specific 
case when W, = Wo. Under such circumstances D sin Wat and F cos w at 
require L' Hospital's rule for proper evaluation, but eventually yield a drift 


term which is a direct function of time. 


The other case of interest concerns harmonics. The series terms asso- 
ciated with Xp require L' Hospital's rule for evaluation when 2W4 = Wo. Under 


such a condition drift terms appear which are direct functions of time. 


An example of drift caused by the existence of harmonics such as when 


2W4 = Wo is shown for the term J sin (2w,t = Wot) ; 


If 2W1> Wo » then (20, — Wo) > 0 (A-55) 


aw, [F - TD (2w, — 9) | 


J sin (2w,t cs Wot) = sin (2w,t - Wot) (A-56) 


2 (2w, — Wo) [1+ 72 (2w, - we)” ] 






Let (2w,t — wot) = (2w, —Wo)t = yt (A-57) 
lim J sinyt = o apply L' Hospital's Rule (A-58) 
y>0 0 
d 
‘dy [aw (F - TDy) sin yt] aw, -7Dsin yt + (F-7Dy)(t cos yt) 








[2y (1 + ray?) } (1 + ry’) + y(27“y) 
y 


(A-59) 


14 


SE « 
— =m © 
(ag pe 96 


mat O° Ea —__ = 


—m ss Em a + em 
ee 


oes 


—— «= «— & 
_— 








tee iret = | RCO) oie (A-60) 
y>0 2 1 2 
If terms for Xn, Xg, etc. were evaluated, drift terms would appear for higher 
harmonics. As the series has an infinite number of terms, it seems apparent 
that any two given input frequencies might create some type of harmonic drift, 
but the magnitude of "'a'' associated with such drift may make the consideration 
of the value a folly when compared to the manufacturing tolerances used when 


constructing gyros. 


The case of negligible elastic restraint shows the possibility of drift terms. 
However, for such a gyroscope, the shortened series representing drift is valid 
only when values of a, b, andc are small compared to unity. Furthermore, 
man can only produce a gyro which is not perfect. Hence terms in the series 
which are of less consequence than the manufacturing tolerances will have only 


a theoretical importance and should not enter practical considerations. 
A.2 Series Solution for Finite Elastic Restraint Between 
the Float and Case About the Output Axis (k # 0) 
The general equation is : 
k= aw, cos wt) Xx+xX+ Tx = bw cos Wot — CWo sin Wot + kb sin Wot 
+ ke cos wot (A-1) 


In accordance with the series established in Section 5. 6, the interest lies with 


the following terms : 


X =X, + Xg + Xp +e ee + Xp y + Xom (A-3) 
where: 
KX, + i + TX, e bwo cos Wot — CWo sin Wot + Kb sin Wot + Kc cos Wot 
(A-61) 
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kX + i + TX. (aw, cos wt) xy (A-62) 


(k-aw, cos wt) x + X+47X = bwo cos Wot — CWo sin Wot + kb sin Wot 
+ Ke COS Wot (A-1) 


Let: 
X = X) + Xp» | (A-64) 


where: 


Kx, + x1 + Tx, = bw cos Wot — CWe sin Wot + kb Sin Wot + kecos Wot 
(A-65) 
Substitute (A-64) and (A-65) into (A-1). 
ee 


e oe ve 
(k-aw, cos wt) (x, + Xo) + (x, + Xo) ai (x, 1: Xp) = bw COS Wot — CWo sin Wot 


+ kbsinwot + ke cos Wot (A-66) 
or 
(k - aw, COs wt) Xo + i aL Te = (aw, cos wt) Xy (A-66a) 
Try a solution for (A-64) of: 
xX, = A Sin Wot + Bcos Wot (A -67) 
7 = A We COS Wot - Bwo Sin Wot (A-67a) 
ey = Aws Sin wot - Bws COS Wot (A-67b) 


Substitute (A-67a) and (A-67b) into (A-65). 
kA sin Wot + kBcos Wot + A Wo cos Wot - Bwo sin Wot —-TA ws sin Wot 
-7B we cos Wot = bwo cos Wot — CWo sin Wot + kbsin Wot + ke cos Wot (A-65) 
Equate coefficients of like terms. 
sinwet: kA —Bws - TAw, = —CWy + kb (A-68) 
2 
2 


COS Wot: kB + AWs - TBO, = bwe + ke (A-69) 
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2 
A(k—T Wo ) = Bwo —CWo + kb (A-68a) 


2 
TBw. —kKB + bw, + ke 
A = 2 2 a T Bwo _ KB oe oe Ke (A-69a) 


Substitute (A-69a) into (A-68a) and multiply through by Wo? 


kTBws kB 4 kbw, + KC - 7” Bus + TRB we - rbws - rke we 





2 Bws e cw’ + kbWs (A-70) 


BL TW aie = 7” wo + T.kws - 05 ] 


ie 3 2 aD 2 
= THW5 —Kk C+ T Kew, -CwWo 





—_—_. = 
akr —-[* ede 1+ a 
oP ~2 





(A-71) 

AK-AT wy + CWo — kb pan kb 
B FS =-- ~ ATWo+ Cc - —— (A-68b) 

2D 9 =) 

Substitute (A-68b) into (A-69a) and solve for A: 

a w5 -k) | ax kb ke 
A = ——“— | — -ATw)+ €-—| +b+ —— (A-69) 

2 a oo eo 


ee 








2 
Al eee, ee 


2 








2 
A= TAK— AT we + TCWy - Tkb - Ak 





(k-aw, cos wt) Xo + Xp + T Xo = (aw, cos wt) Xy 


Let: = 
Xp = Xg+ Xq 


where: 


KX, + ss i re = (aw, cos wt) Xy 

Substitute (A-73) and (A-74) into (A-66a). 

(k —- aw, COs Wt) (x. + X4) 4 (x, + x4) + T Ger 4 x4) = (aw, cos wt) Xy 

or 
(k - aw, Cos wt) X, + X, + TX, = (aw, cos wt) Xe 

A solution is desired for the following equation: 

KX + re Tm = aw, Acos wt sin Wot + AW, Bcos wyt cos Wot 
Let: 

XQ = D sin w,t sin Wot + E sin w,t cos Wot 


+ Fcos Wyt sin Wot + Gcos w,t cos Wot 
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+b+ 


ke 
=P) 


(A-72) 


(A-66a) 


(A-73) 


(A-74) 


(A-66a) 


(A-66a) 


(A-74a) 


(A-75) 





Xg = Dw, cos wy4t sin Wot * Dwo sin wt COS Wot 
> Ew, cos wt cos Wot - EWo sin wt sin Wot 
~ Fw, sin w,t sin Wot + Fwo cos w,tcos Wot 
- Gw, sinw,tcos Wot - Gwo cosw,t sin wot (A-75a) 
Xg = -~Dw% sin wt sin Wot + Dw, Wo cos w,t COS Wot 
+DW Wo COS WE COS Wot- Dw5 Sin wt sin wot 
-~Ew? sin wt cos Wot - EW) Wo cos w,t Sin Wot 
~EW, Wo cos w,t sin Wot — Ews Sin w,tcos wot 
-~Fws cos wt sin Wot - Fw Wo sin wt cos Wot 
— Fw, Wo sin wt cos Wot - Fw cos Wy4t sin Wot 
~Gw? COS W,tCOS Wot + Gw Wo Sin w,t sin wot 
+ GW, Wo sin w,t sin Wot ~ Gus cos Wt COS Wot (A-75b) 


Now substitute (A-75), (A-75a) and (A-75b) into (A-74a) and solve for D, E, 
F, and G. 


(kD sin w,t sin Wot + kE sin w,t cos Wot + kKF cos wt sin Wot + kGcos w,t cos Wot) 
+ (Dw, cos w,t sin Wot + Dwo sin w,tcos Wot + Ew, cos Wt COS Wot — EWo sin w,tsin Wot 

-Fu, sin w,t sin Wot + Fwo cos Wy4t cos Wot - Gw, sin W,t cos Wot GW» cos w,t sin Wot) 
+ (-7TD we sin wt sin Wot + TDW Wo cos w,t cos Wot - TDW1%o cos Wt cos Wot 

-TD ws sin wt sin Wot - TE ws sin w,t cos Wot - TEW) Wo cos w,t sin Wot 

—TEW1 Wo cos wt sin Wot ~ TE ws sin Wy4t cos Wot ~TF we cos Wt sin Wot 

-TFW,Wo sin w,t cos Wot _ T FW) Wo sin w,t cos Wot -TF ws cos w,t sin Wot 

—FG we cos w,t cos Wot - TGW1Wo sin wt Sin Wot + TGW Wo sin w,t sin Wot 

-TG ws cos w,tcos Wot) = Aaw, cos wt Sin Wot + Baw, cos wt COS Wot (A-76) 
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Now equate like terms of the various trigonometric function combinations 


and solve for D, E, F, and G. (Note: A and B are known) 
(sin wit sin Wot) 
kD - Ewo ~ Fw, -7TD wy - TD ws + TGW Wo + TGW Wo =a 
kD -Ewo - Fo, - D(r ws + F ws) +27 GW Wo = 0 
D(k- 7 we - 7 we) - Ew - Fw, +G(27w, 0) =,0 (A-77) 
(sin wyt cos Wot) 
KE + Dwo - Gw, —-TE ws -TE ws - TFW Wo - TFW 1 Wo = 0 
DW + E(k-7 wo -T we) = F(27 wo) — Gw, =0 (A-78) 
(cos w,t sin Wot) 
kF + Dw, -Gws -—T EW Wo ~ TEW {Wo -TF wf -TF ws = Aaw, 
Dw, - E(2 T WW) + F(k-T wt -T ws) - Gw, = Aaw, (A-79) 
(cos wit cos Wot) 


kG + Ew, + Fwo + TDW Wo + TDW - TG ws -TG ws = Baw, 


D(27 1 Wo) + Ew, ib Fwo + G(k- 7 ws -T ws) = Baw, (A-80) 
D(k - 7 ws -T ws) = fo G(2 T W4Wo) =0 (A-77) 
Dwy + E(k- Tw - 7 W9) - F(27W1W9) - Gw, = 0 (A-78) 
Dw, + E(2 T W4Wo) ka ws —T ws) — Gw = Aaw, (A-79) 
D(2 715) + Ew, + Fw, + G(k-7 wi) = Baw, (A-80) 


Solve for D, E, F, andG by determinants. 
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0 —Wo —W, 
0 (k - Tw — Tws) -2T W1Wo 
Aaw 2T W4W (k - Tw" - Tw?) 
1 eZ il 2 
Baw, Wy Wo 
i= 
2 2 
Ww (ke tase) —~2TWqW 
2 il 2 2 
w 2TW4W (ke — Tee — Tae) 
1 1° 2 ih 2 
2TW Wo Wy Wo 


2 


(k- Tw! — TH) 


(A-82) 


Let the denominator of the determinant be known as A or (A-82). 


(k — Tw - 


A = (k- Tw — Tw) 2TW1 Wo 


aT 


“ 
+ Wo Wy 


2 TW Wo 


=o 
1 


2 TW Wo 


“2 
7 TW1Wo Wy 


2 T Wo 


2 
Two) —2TwW1wW9 


(k- Tw? - Tw) 


29 


—2TW1Wo 
(k- Tw? - TwS) 


oo 


(k — Tw = Tw) 
2TW1Wo 
aT 
2 2 
(k= TW - T Wo) 
2TW Wo 


a 
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(k= Tw = TW) 


74 


(k= Tw = Twa) 


(k - Tw = 75) 


= 2 W1Wo 
2 2 
i TW - T Wo) 
= 


(A-82) 





o 
A = (k- Twit ~ Tw) | (x 702 703) + ws (k - Tw - ws) -2 704 w5 
+ ws (k- Tw -TwS) + 2 rus we e 4722 ok (kro? - 702) | 
2 
+ Wo | wg(k - 709 70) + ws e 27 we We (k- 7 wy - TW) 


4: 47 * wwe - w Wy + 27 we w, (k-Tw4 — rae] 
2 
2 
=D Tw, ws (k- Tw - 7 we) - ws 4 arated | 
i (-27w, Wo) | 270,08 = w Wolk - 7 ws - Tw) - w Walk - 7 we - 7 wW5) 


2 
+ 27 W4 Wo (k - rw r we) at WW + arsat ud | (A-82) 


4 2 
A = (k-7 w? — 7 WS) is (2% + 2w2) (k- 7 we 7 we) 


+ 87 ws we (k- Tw 7 ws) = 2wit we s ws + w5 —- 16 tw Wo (A-82a) 


Determinant Form of Coefficients D, E, F, and G. 





0 —Wo —Wy 2 TW1Wo 
0 (k- Tw - Tr ws) -27W,Wo — Ws 
Aaw, 2 TW Wo (gaa wi —T ws) —W9 
Baw, Wy Wo (k-7 we -—T we) 
& (A-81) 
(k-7 wy - 7 we) 0 —W4 27 W1Wo 
Wo 0 -2 TW1Wo —W 
Wy Aaw, (k-7T we 7 ws) —Wo 
2TWsWo Baw, Wo (k-—T wet ws) 
E = 
. (A-83) 
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(k — T wy 7 w5) —Wo 0 2TWsWo 


a 2TWWo Wy Baw, (k — rw - TW) 
a (A-84) 
G= 





a (A-85) 


AD = Aaw [2 rw 05 (k- Tw} -7 ws) a 47 wwe + w? ~ 0105 


2 
2 27 wwe (k- Tw —7 w5) a w, (k= 7 wt - Tw) ‘| 
~Baw, [-27 0405 + 27 We wy PGE W1 Wo (K-T ws - Tw) 
3 2 
a 8 Tw ws 409 (k- 7 OF = 7 wR) — wy Wo (k- T wy - 7 09) 


(A-81) 


AD = Aaw, | #79 (k- 7 wt =T we) + 4ar7 ws 


we + w - w105 
(k 2 25° Ba 2TW Ko aa 
2 
= 27 WW (K-7 wt - TH5) ~ 87 wwe + 2w,w9(k- 7 WH ~T Ws) 
(A-81) 
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AE = ~Aaw, | - 270,09 (k- 7 ot 7 w9) + 27 Wi Wy + 27 Ww 


AE 


AF 


AF 


AG 


Se sa3 2 2 2 2 
+8T W1 Wo + W 1 Wo(k - TW; - T Wo ) + W 1 Wo(k — TW, - T Wo) 
+ Baw, | 2 TW (k= rws _ rw) 2s ws ee, TW (k= rw = rw) 
2 
2 
+47 ws we + w(K - rw — TwW5) “ 2103 | (A-83) 


3). 3 


2 
Aaw, | 270, w5(k- rit — rw) -2 TW Wy - 2Tw4ws - 8T W1 Wo 


_ 2W Wo(k - rw - rw5)| + Baw, [ + Tw10% (k - rw - rw) + ws 


2 
2 
+47 we we + w4(k- Tw4 - TS) + 01% | (A-83) 
Aa (k we wy 2 2 ye 2 2 2 
— 47 wh wh (k- Tay Tap) + w} (K-Twy - Twh) + w(K Tw} - TwH)| 


2 
- Baw, [ ~walk- ray - 7409) + ws Wo a 477 wwe 


= 27w w(k-T we -— Tw) + 2 Tw? Wolk - Tw? - TWH) - 05 
(A-84) 


2 
2)" + 47 ws Wo + (4-477 wt Ww 4 ws) 


2 


Aaw, | (k- ret -Tw 
(k ere) + Baw, | w eK —- WW 
~ 1 2 1 2 1 2 2 


2 
477 w ws t wo | (A-84) 


2 
- Aaw, | wo (k - rw _ Twe) + 4 72 wo we - WW 


+ 2 Tw Wo (k — Tw = Tw) 2 Tw Wo (Kk — Tw _ rw) + 03 | 
5 
+ Baw, [ (x - rw - TW) eZ Tw Ws = 2 Tw we + w8 (k - Tw Tws) 
+4 r we ws (k - rw = Tw) + ws (k - rw = 703) | ( A-85) 
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2 
AG = Aaw, | ww, - wa(k- Tw Tw) - 47? wt ws 


3 
- 470% Wo (K — Tw - 7 Ws) - ws + Baw, (k - rw - 7 ws) 
2 OW. 2 2 
+ (w, +47 oP we + w2)(k-7 we -702)] (A-85) 
(k-aw, cosw,t) x, + X, + TX, = (aw, cos w,t) x, (A-66a) 
Let X4 = Xp 4 Xe (A-86) 
Xg = D sin wt sin Wot + E sin w,t cos Wot 

+ F cos w,t sin Wot + G cos wt cos wot (A-75) 
where: KX, + Xe + T Xp = (aw, cos wt) Xe (A-63) 


Substitute (A-86) and (A-75) into (A-66a). 
(k-aw, cos wt) (x, + Xe) + (x, + Xp) + T (X, + Xe) 
= Daw, sin wt cos w,t sin Wot + Eaw, sin Wt cos wyt cos Wot 
“+ Faw, cos” wyt sin Wot + Gaw, cos” w,t cos Wot (A-87) 


Try a solution for Xp where: 


s a Daw, | Eaw 
kX, + Xp + TX, = sin 2w,t sin Wot + 





Agcpes 
sin 2w,t cos Wot 

















Faw, Faw 
+ sin Wot + cos 2w,t sin Wot 
Gaw Gaw, 
7 cos Wot + cos 2w,t COS Wot (A-88) 
2 
Let 
x 


7 Hy sin 2w,t sin Wot + I sin 2w,t cos Wot 
+ J cos 2w4t sin Wot + Locos 2w,t cos Wot 


+ Msin wot + N cos wot (A-89) 
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Xp could be written in an equivalent form with certain terms containing 
(2w,t - Wot) by use of trigonometry relations. Thus a time dependent 
drift term is seen to appear when the harmonic situation 2W4 = Wo 


exists. 


X = 2H, w,cos2w,tsinwot + H, w, sin2w,tcos wot 
+ al W, Cos 2w,t cos Wot -~ I Wo sin 2w4t sin Wot 
- 25 Wy sin 2u,t sin Wot + JW, Cos 2w,t cos Wot 
- 2L Wy sin Zw ,t cos Wot - L Wo COS 2w,t sin Wot 


+ Mw, coswot - Nwo sin wot (A-89a) 


Combining terms of X,, one obtains: 
Xp = (2H), _ Lwo) cos 2w,t sin Wot + (HW. —_ 2Lw,) sin 2w,tcos wot 
(2Iw, + Two) cos 2w,tcos Wot * (Iwe a 23w4) sin 2w,t sin Wot 


+ 


4 


(Mw) COS Wot - (Nwo) sin Wot (A-89a) 


Xe = (2H, = Lwo) (—2w,) sin 2w,t sin Wot + (2H, _ Lw.) (wo) cos 2w,t cos Wot 
i (Ho = 2Lw,) (20) cos 2w,t cos Wot + (H, > = 2Lw ,)(-wo) sin 2w,t sin Wot 
+ (2Iw, +: Jwo)(-2w,) sin 2w,t cos Wot + (21w, e JWo)(-Wo) cos 2w,t sin Wot 


= (Iw + 23w,)(2w,) cos 2w,t sin Wot = (Iw. fe 2Iw )(Wo) sin 2w,t COS Wot 


(M ws) sin Wot —(N we) cos Wot (A-89b) 


Combining terms of Xp, one obtains: 


Xe = (2Lw, Wy - 4H, 04 -H, w5 + 2Lw,w») sin 2w,t sin wot 
+ (2H, 01% -L ws 4 2H W1Wo - 4L ws) cos 2w,t COS Wot 
+ (-4I ws — 2Jw Ws - lws = 2IW 1 Wo) sin 2w,t cos wot 
+ (- 21wW Wo - Jw% _ 21w We - 43 w,) cos 2w,t sin wot 


+ (-Mws) sinwot + (-Nw) COS Wot (A-89b) 
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Now substitute (A-89), (A-89a), and (A-89b) into (A-88) and solve for 
coefficients HyI, J, L, M, and N., 


[ KH, sin 2w,t sin Wot + kI sin 2a) yt cos Wot + kJ cos 2w4t sin Wot 
+ kLcos 2w,t cos Wot + KM sin Wot + KNcos wat| + 

(2H 0, - Lv) cos 2w,t sin wyt + (H, Ww» - 2Lw,) sin2w,t cos wot 
+ (21@, + JW) COS 2w,tcoswot — (I Wo + 23 w,) sin 2w,tsin wot 
+ (M Wo) cos Wot — (N Wo) sin wot | ot 

[7 (2L we - 4H, wy ~ Hy ws + 2Lw, 0) sin2w,t sin wet 

+ 7 (2H, W4 Wo - Lws + 2H, WW. - 4L, wf ) cos 2w,t cos Wot 

+ T( -41w% — 23 04 Wo - Lws - 23 W4 Wo) Sin 2w,tcos wot 

+7( —21W Wo ais ~ 21W, Wo ~ 43 w%) cos 2W,tsin wot 


+7 (~Mws) sin Wot + r(=No¥) cos wot | = 


Daw, Baw, 
sin 2w,4t sin Wot : sin 2w,t cos Wot 
2 




















2 
Faw, Faw, 
+ Sin wot + cos 2w,t sin wot 
2 1 2 
2 2 
Gaw, Gaw, 
+ COS Wot + cos 2w,t cos Wot (A-90) 
2 2 


Now equate like terms of the various trigonometric functions and solve 
for H,, 1, J, L, M, and N. (Note: A, B, C, D, E, and F are known) 


sin 2w,4t sin Wot 





9 9 Daw, 





(A-91) 
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sin 2w,t cos Wot 
Eaw, 
2 





2 2 
kI + (H, %, - 2Lw,) + T(-41@, — 43 w, Wo -1w95) = 


9 9 Baw, 
H, (wy) + (k-47w; -Tw) + J(-47 ww.) + L(-2w,) = 





(A-92) 
cos 2w,4t sin Wot 


2 


9 Faw, 
kJ + (2Hw, - LW») + T(- 41, W -JI wo - 43 w5 


2 


) = 





Faw, 





H (2w,) : I(-47w, wo) +I(k-Tws — 4; we) + L (-w,) = 
(A-93) 


cos 2w,t cos Wot 




















9 9 Gaw, 
foley (2Iw, + Two) a T (2H WW. -Lw 2H WW, - 4L.w5) = 
9 9 Gaw, 
H (47 w,wo) + 1(2w,) + J(wo) + L(k-T ws - 475) = 
Zz (A-94) 
cos Wot 
9 Gaw, 9 Gaw, 
Mw, —TNwy = ; M (wo) + N(-Tw5) = 
(A-95) 
sin Wot 
9 Faw, 9 Faw, 
— Nwo -T MW = 9 M (- Tw5) + N(-wo) = 
(A-96) 
Solve (A-95) and (A-96) for M and N. 
Gaw, 
M = TNw, + (A-95a) 
2Wo 
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Substitute (A-95a) into (A-96). 











~ Nwo — TWo TNW, + = 
2 3 GTAW, Wo Faw, 
= NO, —- Ni}. 2 
2 2 9 


Faw, + TGAW 1 Wo 


2 3 


aw (F a T Go) 


—2Wo(1 + ae ws) 


: ~aw,(F + T GWo) 

















(A-97) 
2Wo(1 f rw) 
Faw, 
N = =% Mw> - (A-96a) 
2Wo 
Substitute (A-96a) into (A-95). 
9 Faw, Gaw, 
Mwy - TW. | -TMw, - = (A-95) 
2Wo 
T Faw,w Gaw 
Midy ¢7 Mus + ee - = : 
2 2 
Gaw, — 7 Faw,w 
M(we, + re oe) = be eee 
2 2 9 
aw (G —T Fwo) 
= ee (A-98) 


2Wo(1 + 7 ws) 
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Daw 1 





m2 2 
(k- 47 w5 - TW) H, + (-wo)1 + (-2w,)J + (47w,W.)L = 























(A-91) 
9 9 aw 
(Wo) Hy + (K-4Twy - Twy)I+ (-47w, Wo) J + (-2w,) L = 1 
4 (A-92) 
2 aw 
(24) H, + (-47w, wo) 1 + (k— Tw ~47wt) J 4 (-wo)L = 1 
e (A-93) 
aw 
(471 Wy) Hy + (20,)1+ (wo) I + (k- Tw ~ 4704) L = 
(A-94) 
Solve for Hi; I, J, and L by determinants. 
Daw, (A-99) 
; (—wo) (~2w,) (4 TW1Wo) 
; (k ~ 4 TW, TWo) C=407 W 1 Wo) (-2w,) 
Gaw 
5 : (2w,) (wo) (k- arw - rws) 
= Z 2 
(k-470, -T Wo) es) (—2,) (47, Wo) 
(Wo) (ke 47W; = rw) (—4 TW Wo) (~2,) 
(2,) (-47W Wo) (k- 470% - Tw) (- w2) 
(4 TW4Wo) (20,) (wo) ic arws = Tws) 


(A-100) 
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Let the denominator of the determinant be called A or (A-100). 









































9 9 Daw, 
= ATW, = T Wo) ; (—2w,) (47 Wo) 
Baw, 
(wo) : (—47 w1 Wo) (—2w,) S 
Faw < 
1 2 2 = 
(2w,) ; (k - 47 w' = T Wo) (—Wo) 
Gaw, 9 9 
(4r W 1 Wo) ; (wo) (k— 47 Wi - T Wo) 
I ne 
r 
9 9 Daw, 
(k _ 47 Ww, - T Wo) (-Wo) 3 (47 wy Wo) 
Eaw 
2 2 
(wo) (k- 47, - T Wo) : | (-2w,) a 
Faw < 
Gaw, 9 9 
(47 w, Wo) (2w,) ; (k-47w, — Wo) 
J = 
Av 
Daw 
(k= 4T ws - Tw?) (—Wo) (—2w,) ; ! 
2 9 Kaw, = 
(wo) (k — ATW - T Wo ) (—47 WW) 5 ° 
(2w,) (—47 w Wo) (k~47w) - Tw) ; 
Gaw 
(47 W4 Wo) (24) (wo) 1 
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(k-4Tw! - TW) (—47 w, wo) (—2w,) 


A= (k-47 w4-T WS) (~47 w, Wo) (k= ar ws - Tw) (—wo) 
(2w,) (ws) (k- 474 - Tw5) 
(wo) (—47 w, wo) (-2w,) 
Ss (-w9) (2w,) (k - 4rw - Tw) (—wo) 
(47, Wo) (wo) (k- 4704 _ Tw) 
(wo) (k-4TwS-TwS) (2) 
+ (—2,) (20) (245 WW) (—Wo) 
(4 Tw, Wo) (24) (k- 47 w% — Tw) 
(wo) (k - ar ws - 7 ws) (-47 W1Wo) 
- (4701 W,) (20) (47,09) (k-4 Tw - TW) 
(4 Tw) Wo) (2w,) (wo) 
(A-100) 


5 
= (k-47 wit - 7 w5) | (k-47 02 - 7 we) + (87 04 ws) + (87 wit ws) 
~(-4w4) (k - 47 w4 - Tw5) -(16 7? w? we) (ke - 47 -T we) 
~ (-w2) (k- 4702 — v8) | 
2 _2,% 2 
+ (Wo) [ (wa) (k - 47.0 — 75) + (1672 ws ws) - (4w) wo) 
- (-8r w4 Wo) (k-47 0% -7w5) - (-8T wtw)(k -47 ws ~ 7 w5) + (w3) | 
e (—2w,) | (-47 w 19) (k -47 ws _ 7 we) A ale 47 ws — 7 we) 
te -8 v3) = (3272 ws WS) = (2w,)(k — 47 we oe — (-20,03)| 
(-4 ) | (-4rw, 08) + (47 0, wok - 47 04 ay". 16 Sw.) 


33 


— (G45 I to,se: 47 wh — Twe) - (20,0 Wo)(k — 47 wt 705) | 


(A-100) 
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4 
— | (k- Twi - Tw) + (Bws - 3277 wt we + 2we)(k-47w? — Twe) 








2 


+ (64 Tw we (K-47 w? T we) + (ws - 1605 + 32 7? we We 





2 2 
+ 1287? wl we + 256wiw, To Bui 5 | (A-100) 
pea (k-4Tw%- Tw2) (47 w,W9) (-2 
aes —4Twy- Two)  (-4T Wy Wo 201) 
2 
(-47 wo) (k — 47 ws = Tw) (—wo) 
(20) (Wo) (k -—47 we —T ws) 
aa (—Wo) (—2w,) (47 w, Wo) 
1 2 2 
; (= 4T W1 Wo) (k-47 04 - Two) (-—wo) 
2 
(2w,) (wo) (k — 47 Wi - rw) 
(-wo) (-2w,) (47 W 1 Wo) 
Ae 5 (k-47w, - Two) (— 47 w, wo) (—2w,) 
(2) Cm) (k- 47 wt — Tw?) 
Gaw (—wo) (—2w)) (47 wo) 
: 1 (k- 479 - Tw) (—47 w, wo) (—2w,) 
2 (—47 WW») (k - 47 ws — we) (—Wo) 


(A-99) 
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Daw 3 
2 
AH, = ( ; 7 | (k- aru? Tw) + (87 Ww; we) + (87 w4 ws) 


2 Z 
+: (45 (k - 47 W5 = 75) = (1672 we ws (Kk _ ar wi _ Tw) 


- (w5) (Kk = 47 w s rap) | 





Eaw 2 
1 2 2 2 2 23 


+ (87 004 Wo) (K — Ar ws _ 7 ws) + (87 w4 Wy) (k- 4rws = 7 we) + (03) | 


Faw 
1 2 2 2 3 
( : Kec ~ T Wo) + (8w)) 





2 
+ (47 w,w5)(k - 47 wi _ Tw) + (3277 ws w5) + (2w,)(k -ATw4 - Two) 


2 (2, ws | 


2 





Z 
+ (647 Fw! w>) ~ (401 w(K - 47 ws - i oD 


Gaw, 4 3 3 2 2 
T W4 Wo) + (167 w} wo) — (47 w, Wo)(k - 47 Wy _ T Wo ) 


(A-99) 
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Daw, 9 D 
| (k - aru? rut)” + (16 TwS w Wo) + 


+ (407 = 1677 we ws - we (k - aT we — “4 + 








Eaw, 9 3 
; (wo)(k - 47 ws - rue) ~ (404 Wo) + (167% w we we ) + 


+ (167 w4 wo)(k - 47 we Tw) + 2 + 


Faw, 9 é 2 eee 
; (20,) (k - 4Tw ‘= 7 we) ~ (2w, ws) + 327 Wy Wo + 





+ (87 W105) (k~ 47w4 - Tws) + (8 2 + 





Gaw, 9 . 2 cae 
; (—47 Ww, Wo) (K-4tw>-Two) + 47W WS + 


+ 647 ws ws Wo + (47 w, Wo) (k - ar ws _ rw) + 


3 
+ 167w, “| 


(A-99) 
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Daw, 
I= = (wo) (—47w, Wo) (~2w,) 











é 2 2 
(20)  (k-4Tw) - Tw5) (—wo) 
(47 W4Wo) (wo) (k —4g we ~ Tw) 
Baw, 9 9 
+ 3 (k - 47 w = TW) (~2w,) (47 w, Wo) 
(2w,) (k-—47 we _ Tw) oe 
(47 w, Wo) (wo) (k — 47 wo, Tw) 
Faw 
; ' (k- 47 ws ~ Tw?) (-2w,) (47 wW Wo) 
(wo) (-47 w, Wo) (-2w,) 
(47 W1Wo) (wo) (k e: Tw ~ Tw) 
Gaw, a 9 
- ; (k -47w; -TW5) (-2w,) (47 Wo) 
(wo) (-47T WW) (-2w,) 
(2w,) (k ~ 47 ws _ 7 w5) (wo) 
(A-101) 
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ew 2 2 23 
vie : (Wo) (k - ‘nn + (167 3) - (40% w Wo) + 


_— 





+ (8rw, Wo) (ee Ar w - 7 ws) + (87 ws Wo)(k-47 w4- 7 ws) = 











3, | 
- eu = 
Eaw 3 
+ 5 t (k- 47 w% - Tw) a. (87 wi we) + (87 we ws) + 
- (1677 w? ws) (k = 4704 - TW) + (409 )(k- 47 - 7) + 
= (ofretra ra) : 
Faw, 9 9 2 3 
- ; (-47 W Wo (kK - 47 wy - TW) + (16Tw) Wo) + 
+ (47 103) = (647° ws ws) + (2w , Wo)(k-47 w4 _ 75) + 
+ (2w,0)(k-47 w* ~ tw")| + 
LZ I 2 
Gaw 
i 2 Z 2 3 
i" 5 (47 WwW) (k- 47 wy = T Wo) - (8w5) “ 


+ (47 105) (k - ar ws ~ 7 W5) + (3277 we ws) = (2w ws) - 
2 
a (20 M(k-47 w4 = 705) | 


(A-101) 
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Daw 











Gaw 





- 


£ 
t 


(wo) (k — aT wt ~ T ws) (-2w,) 
(2,) (=47 W4 Wo) (—Wwo) 
(47 W4Wo) (2w,) (k 247 ws - Twe) 
(k ~47 ws ~ Tw) (~Wo) (47 W4Wo) 
(2w,) (—47 w Wo) (—Wo) 
(47 w, Wo) (20) (k- 47 wf - Tw?) 
(k — 47 ws _ Tw) (wo) (47 W 1 Wo) 
(wo) (k ~4T we ~ T we) (—2w,) 
(47 w, wy) (2) (k- 47 w2 - Tw) 
(k-47 ws Tw) (—Wo) (47 w wo) 
(wo) (k ~ 47 wi _ Tw) (—2w,) 
(2w,) (~47 w, Wo) (-Wo) 
(A-102) 
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Eaw 








"| 





Faw 
+ 
2 


{ 


Gaw 


2 


j 





(Wo) 
(2w,) 


(47 w Wo) 


2 2 
(k= 47 ws - T Wo ) 
(2) 
(47 w1 Wo) 


(k— 47 wi ~ 7 ws) 
(wo) 
(47 W1Wo) 


(k- 47 ws _ 7 w5) 
(wo) 
(20,) 


(k — Tw 7 Tw) (-47 w wo) 


(-47 w Wo) (k 47 w4 - 7 w5) 


(—W5) (-2w,) 
(~47 wo) (k- 47w4 - Tw) 

(—wo) (-2w,) 

(k-47 Ww - 7 WS) (—47 W, Wo) 
(20,) (wo) 
(—Wo) (-2w,) 

(k ~ 47 we -- rw5) (-47 W 1 Wo) 

(-47 w, Wo) (k- 47 w% - Tws) 
(A-103) 
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XJ and AL are left in the form of a third order determinant. Fur- 
ther determinant breakdown will yield terms similar to those associated 


with AH and Al, 


The determinant form does permit the detection of the common 
factor "a" in all the coefficients. As in the case for zero elastic re- 
straint, terms beyond X_ may be neglected if the order of magnitude 
of ''a" is much smaller than unity. Each successive term in the series 
of x is smaller than its predecessor in accordance with the value of 
bal. 


K=X,+XgtXe+... . + Xp 4 + Xo (A-3) 
x, =A sin Wot + Bcos Wot (A-67) 


Xg =D sin wt sin Wot + E sin wt cos Wot 


+ F cos wt sin Wot + G cos Wyt cos Wot (A-75) 


Xp = H, sin Zw t sin Wot + I sin 2w,t cos Wot 
+ J cos 2w,t sin Wot + Locos 2w,t cos Wot 


+M sin wot + N cos wot (A-89) 


Drift terms may be shown to exist by applying trigonometry re- 
lations to Xa and Xp Equal input frequencies and harmonics will 
cause time independent terms of constant angular displacement to 
appear but there are no time dependent terms as in the case of zero 


elastic restraint. 
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Xq = . [ cos (wt — wot) - cos (wt es Wot) | 
+ a [ Sin(w t+ Wot) + sin (wt — Wot)] 
- _ [ sin (wt + Wot) _ sin (w,t - Wot) | 
ue 2 [ cos(wyt + Wot) + COs (w,t ns Wot) | (A-75a) 
Xp = = [ cos (2w,t — Wot) - cos (2w,t + wot) | 
' e [ sin(2w,t + Wot) + sin (2w,t & Wot) ] 
r + [ sin(2w,t + wot) — sin (2w,t-wot) | 
" = [ cos (2w,t + Wot) + COs (2w,t = Wot) | 
+ Msin wot + N cos wot (A-89a) 


Xy terms are pure Sinusoids and the time average aang at the end 
of complete cycles is zero. Xa terms are sinusoids which do not 
display time dependent characteristics if Wy = Wo. Constant dis- 
placements such as D/2 and G/2 do appear when W1 = We. Likewise 
in X5 terms, time dependent terms are lacking when 204 = Wo but 
constant angular displacements such as H/2 and L/2 are present. 
This phenomena can be easily explained through the coefficients 
attached to the sinusoids. The denominators of such coefficients 
become zero for W1 = Wo and 2W4 = Wo in the case for zero elastic 
restraint. However, a finite elastic restraint keeps a real value 
attached to the coefficients and thus indeterminate forms do not ap- 


pear when sin(w,t = Wot) and sin (2w,t — Wot) are zero. 
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This analysis is verified by Ref. 10 for the case where Wy = Wo 


and k is finite. 


One other possibility is seen to exist. It concerns A and ), 
the denominators associated with the determinant forms of the co- 
efficients. If A and A can be made equal to zero by some combination 
of W4 and Wo; drifts could possibly occur. Attempts to solve A for 
Wy to give A = 0, lead to “an equation in Ws of the eighth order. Simi- 
lar difficulty appears whet » is set equal to zero and an attempt is 
made to solve for Wy: With such values of Ws and Wo which satisfy 
A=0 ord =0, L' Hospital's Rule could be applied and the possibility 
of a drift rate could be investigated. This appears as a field for more 


complete investigation. 
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